Apparatus and method using measurements taken while drilling to generate and map mechanical boundaries and mechanical rock properties along a borehole

ABSTRACT

The present disclosure involves a novel way of using drilling vibrations generated by the deformation of a rock formation in response to forces acting on the rock formation, where the forces are related to a drill bit and/or drilling fluid system, to identify the nature and occurrence of fractures, fracture swarms and other mechanical discontinuities (boundaries) such as bedding planes and/or faults that offset or otherwise separate rock formations with different mechanical rock properties.

RELATED APPLICATIONS

This application claims priority under 35 U.S.C. §119 to U.S. Provisional Patent Application No. 62/216,786 entitled “APPARATUS AND METHOD USING MEASUREMENTS TAKEN WHILE DRILLING TO MAP MECHANICAL BOUNDARIES AND MECHANICAL ROCK PROPERTIES ALONG A BOREHOLE,” filed on Sep. 10, 2016, which is incorporated by reference in its entirety herein for all purposed. This application is also a continuation-in-part of U.S. Nonprovisional patent application Ser. No. 15/182,012 entitled “APPARATUS AND METHOD USING MEASUREMENTS TAKEN WHILE DRILLING TO MAP MECHANICAL BOUNDARIES AND MECHANICAL ROCK PROPERTIES ALONG A BOREHOLE,” filed Jun. 14, 2016, which is a continuation of U.S. Nonprovisional patent application Ser. No. 14/850,710 entitled “APPARATUS AND METHOD USING MEASUREMENTS TAKEN WHILE DRILLING TO MAP MECHANICAL BOUNDARIES AND MECHANICAL ROCK PROPERTIES ALONG A BOREHOLE,” filed Sep. 10, 2015, which claims priority under 35 U.S.C. §119 to U.S. Provisional Patent Application No. 62/048,669 entitled “APPARATUS AND METHOD USING MEASUREMENTS TAKEN WHILE DRILLING TO MAP MECHANICAL BOUNDARIES AND MECHANICAL ROCK PROPERTIES ALONG A BOREHOLE,” filed on Sep. 10, 2014, all of which are hereby incorporated by reference in their entirety for all purposes.

TECHNICAL FIELD

The present disclosure involves measurement while drilling techniques that provide mechanical rock properties, and from which fractures and other mechanical boundaries may be identified and used to improve drilling and completion practices, among other things.

BACKGROUND

The production of hydrocarbons (oil or gas) can be generally distilled into two primary steps—drilling a borehole to intersect hydrocarbon bearing formations or oil and gas reservoirs in the subsurface, and then completing the well in order to flow the hydrocarbons back to the surface. The ability of a well to flow hydrocarbons that are commercially significant requires that the borehole be connected to oil and gas bearing formations with sufficient permeability to support the flow rates that are needed to account for the costs of developing the field. In many instances, however, commercially viable flow rates cannot be obtained without the use of various advancements including horizontal drilling and hydraulic stimulation due to the type of formation or reservoir being developed.

More specifically, unconventional resource plays are areas where significant volumes of hydrocarbons are held in reservoirs with low primary permeability (nanodarcy to microdarcy) and low primary porosity (2-15%) such as shales, chalks, marls, and cemented sandstones that generally do not have sufficient primary permeability to yield commercial quantities of hydrocarbons. Compared to “conventional” reservoirs, unconventional reservoirs have a much lower hydrocarbon density per unit volume of rock and much lower unstimulated hydrocarbon flow rates, making commercial development impossible without hydraulic stimulation of the reservoir rock fabric. Fortunately, unconventional reservoirs are often regionally extensive covering thousands of square miles and containing billions of barrel of oil equivalent (BOE) of potentially recoverable hydrocarbons.

The economically viable production from unconventional resources has only been made possible by the improvement and combination of horizontal drilling, wellbore isolation, and hydraulic fracture stimulation treatment technologies, among other techniques. Generally speaking, horizontal drilling involves first vertically drilling down close to the top of the unconventional reservoir and then using directional drilling tools to change the orientation of the wellbore from vertical to horizontal in order to contact greater areas of the reservoir per well. The term “horizontal” drilling as used herein is meant to refer to any form of directional (non-vertical) drilling. Horizontal drilling, although having been performed for many decades prior to intensive unconventional resource development in the early 2000's, has been evolved to provide cost effective provisioning of the long horizontal borehole sections (5,000′ to 10,000+′) required to contact commercially viable volumes of hydrocarbon bearing reservoir rock. Hydraulic fracture stimulation involves pumping high volumes of pressurized fluid into the borehole and through targeted perforations in the wellbore to create large networks of cracks (fractures) in the formation that create enhanced reservoir permeability and so stimulate greater quantities of oil and gas production. Proppant is usually pumped along with the fluid to fill the fractures so permeability is maintained after the pumping is stopped and the fractures close due to reservoir stresses. Proppant can range from simple quarried sand to engineered man-made materials.

Isolation generally involves the use of some form to technology to focus where fracturing occurs at specific locations along the well bore rather than stimulating the entire length of an open wellbore. In the development of unconventional resources it is desirable to drill horizontal wells perpendicular to the direction of maximum horizontal compressive stress, because hydraulically induced fractures will grow primarily in the direction of maximum horizontal stress. When the wellbore is oriented perpendicular to the maximum horizontal compressive stress, this geometry allows for the shortest, and hence least expensive, well bore length for the volume of reservoir stimulated.

Rapidly evolving wellbore isolations techniques, such as swellable packers, sliding sleeves, and perforation cluster diversion have all assisted in reducing the cost of isolating sections of the wellbore for more targeted and more concentrated hydraulic stimulation. Hydraulic fracture stimulation has been utilized on low permeability wells for decades as well. But the use of low viscosity, simple fluids pumped in very high volumes and rates, and with large volumes of associated proppant, has been the most important aspect of contacting the greatest amount of low permeability, low hydrocarbon density reservoir rock.

Various suites of drill string or wireline conveyed well logs such as dipole sonic or natural fracture image logs can identify and quantify this variability on a scale that is useful to completions design, but existing tools are currently too expensive to run on anything but a very small fraction of unconventional wells drilled. Conventional techniques, such as dipole sonic and natural fracture image logs, are based on inferred information and not involved directly measuring the interaction of the drill bit with the formation. Instead, dipole sonic involves the transmission of acoustic signals (waves) from a controlled active acoustic source, through the rock formation in the areas of the well bore to a receiver typically several feet from the source, to measure the velocity of the waves through the formation. Natural fracture image logs involve measuring the resistivity of the formation along the walls of the wellbore. Natural fracture logs are of limited use in wells using oil based mud, which has an inherently high resistivity and masks some fractures. These techniques are often cost prohibitive and limited in effectiveness. As a result, almost all wells are completed using geometrically equal spacing of zones isolated (referred to as stages) and stimulated. Thus, for example, hydraulic fracturing is inadvertently performed routinely on individual stages with significantly varying rock properties along the isolated section, resulting in the failure to initiate induced fractures in less conducive rock and so potentially bypassing substantial volumes of hydrocarbon bearing rock. In such instances, post stimulation testing of individual zones or stages shows that a significant percentage of the hydraulically stimulated zones are not contributing to hydrocarbon production from the well. Variations in the density, size and orientation of natural fractures can have a major influence on overall well initial production, long term decline rates, and stage to stage contributions. Formation hydrocarbons are transported from the rock matrix to the producing wellbore through some combination of induced hydraulic fractures and natural occurring in-situ fractures.

Currently, less than 1% of all wells drilled and completed have suitable data to adequately quantify reservoir heterogeneity on a scale that can be used for targeting individual stimulation intervals.

It is with these observations in mind, among others, that aspects of the present disclosure have been conceived and developed.

SUMMARY

Aspects of the present disclosure involve a method of characterizing rock properties while drilling comprising: receiving acoustical signals obtained from one or more sensors positioned on a component of a bottom hole assembly where the sensors (e.g., accelerometers or strain gauges) are in operable communication with at least one data memory to store the acoustical signals. The acoustical signals, which may also be considered vibrations, are generated from a drill bit interacting with a rock formation while drilling a wellbore. The method further involves processing the acoustical signals to obtain at least one set of data values representative of a mechanical rock property of the rock formation along the wellbore created by the drill bit interacting with the rock formation for a period of time. The method may further involve aligning the data values, or mechanical rock properties generated therefrom, with an axis of material symmetry. With the aligned data, for example, the method may identify elastic coefficients of C₅₅ as well as C₁₁, C₁₂, C₁₃ and C₃₃. The method may further involve identifying a change in the at least one set of data values where the change representative of the drill bit crossing a mechanical rock property discontinuity while drilling the wellbore.

Another aspect of the present disclosure involves a method of characterizing rock properties comprising obtaining, using a processor, dynamic signals from one or more sensors positioned on a component of a bottom hole assembly, the dynamic signals generated from a drill bit interacting with a rock formation while drilling a wellbore. The method further involves aligning, using the processor, the dynamic signals in relation to a position of the drill bit with respect to an axis of symmetry of the rock formation and processing, using the processor, the aligned dynamic signals to obtain at least one set of data values representative of a mechanical rock property of the rock formation along the wellbore created by the drill bit interacting with the rock formation for a period of time.

Yet another aspect involves a method of characterizing and mapping rock properties along a well bore including receiving acoustical signals obtained from one or more sensors positioned on a component of a bottom hole assembly, the acoustical signals generated from a drill bit interacting with a rock formation while drilling a wellbore, aligning the acoustical signals in relation to a position of the drill bit with respect to an axis of symmetry of the rock formation, and processing the aligned acoustical signals to obtain forces acting on the drill bit interacting with the rock formation while drilling the wellbore and to obtain displacements of the drill bit interacting with the rock formation while drilling the wellbore. The method further involves processing the forces and the displacements to obtain at least one set of data values representative of a mechanical rock property of the rock formation along the wellbore created by the drill bit interacting with the rock formation for a period of time.

Another aspect involves an apparatus comprising a processing unit in communication with at least one tangible machine readable media including computer executable instructions to perform the operations of the various methods, or portions thereof, discussed herein. The processing unit may be part of a measurement while drilling apparatus and/or a computing device configured to process data collected from the measurement while drilling apparatus.

These and other aspects are disclosed in further detail in the description set out below.

BRIEF DESCRIPTION OF THE FIGURES

FIGS. 1A-1B illustrate reservoir-to-well connectivity where brittle rocks are generally associated with larger fracture creation and better proppant support that is more permeable than ductile rock that produces smaller, less productive fractures which are prone to rapid compaction and closure and are less permeable.

FIG. 2A-2B is a diagram of a drill bit assembly including sensors for measuring bit accelerations and forces on the bit, and which includes at least one processing unit and tangible storage media in which to store acceleration and/or force data, and which may also store processed acceleration and/or force data of the drill bits interaction with a formation while drilling.

FIGS. 3A-3C illustrate the cutting action of a drill bit, stick slip fracturing of a rock formation, torque-on-bit and weigh on bit forces, and related torque and displacement curves.

FIG. 4 illustrates root mean square computations for axial acceleration of the bit as used in computing rock properties.

FIG. 5A is a rotary displacement spectra obtained from measurement while drilling data;

FIG. 5B is a diagram illustrating lateral and rotary acceleration of the bit, while cutting, and useful in computing mechanical rock properties among other advantages;

FIG. 6 is a diagram illustrating drill bit behavior when encountering a mechanical discontinuity or geological boundary, the cutting face of the drill bit may change its orientation, as detectable from acceleration data, in response to the orientation and stresses acting on the heterogeneity.

FIG. 7 is a graph of rock strength curves computed from torque and penetration per revolution information based on measurement while drilling data.

FIG. 8A illustrates stress strain relationships based on an orientation of a well relative to a transverse isotropic axis of material symmetry and mechanisms whereby measurement while drilling techniques may be used to compute elastic coefficients, from which temporal and/or spatial variations in one or a combination of more of the measurements obtained from the geophysical signal processing techniques are used to identify the nature and occurrence of fractures, fracture swarms and other mechanical discontinuities (boundaries) such as bedding planes and/or faults that offset or otherwise separate rock formations with different mechanical rock properties.

FIG. 8B is a stress strain curve reflecting the relationships of FIG. 8A.

FIG. 9A illustrates a constitutive stress strain relationship for an axis of material symmetry parallel a well axis, and computations using torque-on-bit, weigh on bit and axial displacement to obtain elastic coefficients from measurement of bit vibration and/or forces acting on the bit.

FIG. 9B illustrates a constitutive stress strain relationship for an axis of material symmetry perpendicular a well axis, and computations using torque-on-bit, weigh on bit and rotary displacement to obtain elastic coefficients from measurement of bit vibration and/or forces acting on the bit.

FIGS. 10A-10C a constitutive stress strain relationships for an axis of material parallel and perpendicular a well axis, using acceleration data.

FIG. 11 illustrates constitutive stress strain relationships for an axis of material symmetry parallel a well axis and axis of symmetry perpendicular a well axis, with Poisson's ratio and Young's Modulus of Elasticity computed from measurement while drilling data of torque-on-bit, weigh on bit, and acceleration data to obtain lateral and axial displacement.

FIG. 12 illustrates stress strain curves obtained by the methods discussed herein.

FIG. 13 illustrates two relative curves of Poisson's ratio for a vertical well, with cross over points identifying likely preexisting fractures in the fracture detection log.

FIG. 14A illustrates two relative curves of Young's Modulus of Elasticity for a horizontal well, and 14B illustrates two relative curves of Poisson's ratio for the same horizontal well, with cross over points identifying mechanical rock properties indicative of preexisting fractures.

FIG. 15 is a diagram illustrating fracture reactivation, if the forces acting on the formation in connection with the drill bit and drilling fluid system when conducting drilling operations are sufficient to overcome the failure criteria of a pre-existing fault, then the fault will slip or fail. Reactivation of a fault or pre-existing fracture can be evidenced by extracting a signal from the drilling vibrations that is related to a microseismic event with attendant primary, compressional (P) and secondary, or shear (S) arrivals. In the case where the fault is perpendicular to the trajectory of the wellbore the P-wave arrival is related the particle motion parallel to the axis of the drill string and the S- or transverse wave is the particle motion parallel to the lateral and torsional motion of the drill string. Deviations of the particle motion of the P- and S-waves relative to the orientation of the well can be used to determine the orientation of the fault.

FIG. 16 is a method of obtaining mechanical rock properties of a formation proximate a well bore from measurements of bit behavior taken while drilling.

FIG. 17 is a method of obtaining mechanical rock properties of a formation proximate a well bore from measurements of bit behavior taken while drilling, where the method involves rotation to align measurements with an axis of symmetry of the formation.

FIG. 18 illustrates a cartesian coordinate system with three angles representing potential rotation angles about each principle axis, where any arbitrary rotation of a compliance tensor into a reference frame such that the principle stresses and strains align with the reference frame can be achieved by a combination of rotations about the three principle axes rotate.

FIG. 19 is a diagram illustrating a drill bit striking a formation at an oblique angle.

FIG. 20 illustrates partial rotations of a drill bit, the partial rotations illustrated by θ₁, θ₂, θ₃, and θ₄.

FIG. 21 is a special purpose computer programmed with instructions to execute methods discussed herein.

DETAILED DESCRIPTION

The present disclosure involves a novel way of using drilling vibrations generated by the deformation of a rock formation in response to forces acting on the rock formation, where the forces are related to a drill bit and drilling fluid system, to identify the nature and occurrence of fractures, fracture swarms and other mechanical discontinuities (boundaries) such as bedding planes and/or faults that offset or otherwise separate rock formations with different mechanical rock properties.

The techniques and measurements from this disclosure are made using downhole tools that are simple and rugged, allowing for a magnitude in order reduction in logging cost to characterize near-wellbore rock mechanical properties and intersected existing fracture locations. The low cost to log a well, typically less 0.5% of the total well cost, allows for widespread use of the technique. Detailed knowledge of rock property variability along a wellbore allows for grouping like-for-like rock types in variable length stages, avoiding losing reserves due to a lack of fracture initiation relative to mixed rock strength stages. Also, knowledge of existing fractures will improve well economics overall as fractures can be targeted for stimulation to improve initial production if appropriate, or avoided for example when setting swell packer locations.

Further elaboration of the method describes how mechanical rock properties, and in particular elastic coefficients of a rock formation, can be determined through the application and use of innovative, new stress-strain relationships that systematically relate measurements of the forces acting on rock formation in connection with the drill bit and drilling fluid system (stress) to the variations in the drilling vibrations generated by the deformation of the rock formation in response to the cutting action of the bit (strain).

The elastic coefficients can be used, in general to describe the deformation of a rock formation in response to the forces acting on the rock formation, and in particular, to predict the deformation of a hydrocarbon bearing formation in response to the forces acting on the formation where the forces are fluid pressures generated during the emplacement of hydraulic fractures in connection with a hydraulic fracture stimulation treatment along a horizontal well.

The processing of drilling vibrations when recorded using sensors deployed in a borehole in connection with a bottom hole assembly (BHA) according to the method disclosed here, can provide measurements of mechanical rock properties including the nature and occurrence of mechanical discontinuities, such as pre-existing fractures, which can be used to target sections of the well where the rock properties are conducive to economical hydraulic stimulation and to avoid sections that are viewed as sub-commercial, where the rock properties are not conducive to economical hydraulic stimulation.

In another embodiment of the method, the elastic coefficients and variations in the elastic coefficients that are obtained while conducting drilling operations can be used for assisting, in real-time, the steering of the bottom hole assembly in order to maintain the tracking of the drill bit through geological formations as are targeted according to the desired mechanical rock properties, especially where the mechanical rock properties are relevant to the production of commercially significant hydrocarbons using hydraulic fracturing stimulation techniques.

The present disclosure generally relates to the production of commercially significant hydrocarbons from oilfield drilling operations and completion operations. Recent unconventional resource development has identified a need for economical determination of rock properties and natural fracture swarm locations along a horizontal well in order to optimize the location and intensity of hydraulic stimulation treatments.

The techniques described in this disclosure will provide new information for selecting hydrocarbon bearing zones by differentiating between brittle rocks generally associated with larger fracture creation and better proppant support that is more permeable than ductile rock that produces smaller, less productive fractures which are prone to rapid compaction and closure and are less permeable. Natural fracture identification also refines the process of hydraulic stimulation optimization by providing direct measurement of zones that offer higher permeability and higher hydrocarbon productivity.

Physical Basis of the Method

Aspects of the present disclosure involve methodologies that use broad band measurements (e.g., continuous, high resolution) of drilling vibrations and drilling dynamics data taken proximate the drill while conducting drilling operations bit to log the mechanical properties of a rock formation.

Drilling vibrations generated by the deformation and failure of a rock formation are generally related to the mechanical properties of the rock being drilled. It is generally understood that the depth of cut or the tooth penetration into the rock is inversely related to the strength of the rock. Higher amplitude drilling vibrations occur in rocks that undergo a greater depth of cut and deeper tooth penetration in response to the forces acting on the formation in connection with the drill bit and drilling fluid system, whereas lower amplitude drilling vibrations occur in rocks that undergo relatively lower depth of cut and lesser tooth penetration. Increased depth of cut indicate the bit is moving into an area of lesser relative mechanical rock strength, and decreased relative drilling vibrations indicate the bit is moving into an area of greater relative rock strength all other things being equal.

Generally speaking, rock formations that take a relatively long time to drill through or where the rate of penetration is slow are generally referred to strong or hard formations and have a lesser depth of cut in relation to rock formations that are relatively weaker and less rigid. These basic principles have enabled the application and use of techniques that take measurements of the hardness of a rock formation by forcing a tool into a rock to make an indentation where the depth of the indentation relative to the force applied is used to obtain a hardness characteristic that is essentially a mechanical property of a rock formation.

The presence of mechanical discontinuities, such as pre-existing fractures, and geological boundaries, such as faults, in a rock formation generally act to weaken the rock formation. Fractured rock formations are generally weaker and less rigid than intact, unfractured or stiff or otherwise competent rock formations. As the drill bit encounters fractures in a rock formation the tooth penetration or depth of cut and subsequently the drilling vibrations will increase, because the rock formation is less rigid because it has been weakened by the presence of fractures. Stated differently, as the drill moves into and through existing fractures, the measured mechanical rock strength will decrease relative to the same rock formation without a fracture or with lesser fractures, for example.

General Description of Fracture Identification in Relation to the Method

Signal processing techniques are used to process the drilling data to identify locations where the changes in the drilling vibration indicate that the drill bit has encountered a mechanical discontinuity or geological boundary. If the changes in the drilling vibration as expressed through the results of the geophysical signal processing techniques are rapid and discrete in both space and time, and then return back to a long-term trend or the levels that were recorded prior to the change in the drilling vibrations, then it would indicate that the drill bit has encountered and crossed a discrete mechanical discontinuity because mechanical rock properties that are discrete in both space and time are uniquely separated from the mechanical properties of a rock formation such as would be in the case of a drill bit penetrating a fracture face. If the changes in drilling vibration are rapid and discrete and continue over a short interval, then that would indicate multiple fractures or a swarm of fractures has been encountered.

If the signal processing techniques indicate that the changes in the drilling vibration are rapid, but then do not revert back to the level prior to the change and instead carry on at a new, significantly different level, then that indicates a mechanical boundary where the mechanical boundary that separates or offsets two different rock formations such as a bedding plane and or fault has been encountered and crossed. Whether or not the boundary is related to a fault or a bedding plane depends on the inclination of the bit with respect to the orientation of the stratigraphy of the rock formation being drilled. Other information may also be used to determine whether or not the mechanical boundary was a bedding-plane fault or bedding plane that acted as a zone of weakness that had experienced measurable displacement in the past.

The description provides a method to evidence the presence of fractures, fracture swarms and other mechanical discontinuities such as faults and bedding planes that offset or otherwise separate rock formations with different rock properties. The approach uses geophysical signal processing techniques that are sensitive to changes in the drilling vibrations where the changes are relative to some baseline, such as a normalized preceding set of drilling vibration data, and whether or not the changes are discrete and then return back to the level prior to the change or are maintained at a new level that is different than the level observed prior to the change.

The following method disclosed below further elaborates on the outlined principles to provide a general, independent method to specify the mechanical properties of a rock formation by processing the drilling vibrations in relation to the forces acting on the formation in connection with a drill bit and, in some instances, a drilling fluid system, which includes the mud motor and the drilling fluid, including mud, that turns the motor. This method specifies the mechanical properties of a rock formation through the application and use of innovative, new stress-strain relationships, among other advances.

Because the cutting depth or the penetration of the bit tooth into the formation is a measurement of displacement, the drilling vibrations describe the strain experienced by a rock formation in response to the cutting action of the bit where greater cutting depths and a greater penetration of the bit tooth relative to the same volume of rock result in a higher strain. Strain is understood to describe a change in volume of a rock under some force. In some aspects, strain is reflected by axial displacement of the bit per turn of the bit. In one example, axial displacement is computed as a double integral of the accelerometer data (axial) for one revolution of the bit to yield a distance measurement for one revolution (or some other known number of turns). If this strain is calculated with respect to time, then the drilling vibrations can be used to describe the strain rate. The converse with regards to the strain and the rate of strain is also held to be evident.

Through the techniques described herein, the drilling vibration characteristics, which may be supplemented with drilling dynamics data including forces on the bit, such as torque-on-bit and weight-on-bit, are translated into mechanical properties of a rock formation. The depth of cut, as obtained based on vibration assessment, may be normalized against direct weight-on-bit and/or torque-on-bit measurements, or weight-on-bit or torque-on-bit measurements extrapolated from vibration information.

General Stress-Strain Relationships

Stress-strain relationships are established by systematically relating forces acting on the formation in connection with the drill bit and drilling fluid system to the geophysical signal processing of drilling vibrations generated by the deformation of the rock in response to the cutting action of the bit. This approach allows elastic coefficients (K) to be derived in accordance with the following equation where (e) is the general deformation (strain) of a rock formation in response to the forces acting on a rock formation (S) (stress).

S=K e

In accordance with some methods set out herein, strain (the motion or displacement of the bit is obtained by signal processing the drilling vibrations as are transmitted along the drilling assembly by sensors deployed in a borehole in connection with a BHA. Stress, in accordance with some methods set out herein, is obtained from either (i) downhole or (ii) surface measurements of torque and/or weight-on-bit, or (iii) in another example, the accelerations of the bit are related to forces on bit where it understood that the acceleration is a representation of force per unit mass. It should be appreciated that forces can be converted to stresses with knowledge of the effective contact area of the bit and formation, and the effective rock volume the bit is acting on. Conversely, forces can be substituted for stresses with the understanding that a geometric correction in relation to the effective contact area is required to obtain absolute values for the mechanical rock properties. One example of such a contact area is the area of the bit.

The equations of linear elasticity are useful for describing the relationship between the changes in shape and position of a material in relation to the forces acting on the material. Such stress-strain relationships are known in general as Hooke's law where the coupling of the stress-strain relationship behavior is described through a matrix of coefficients whose values depend on the conditions used to load the material in relation to the structural symmetry of the material being loaded. These coefficients are colloquially known as the cij's and can be arranged in well-known and convenient forms to represent Young's Modulus of Elasticity (YME) and Poisson's Ratio (PR). In one example of the technique presented here, the YME and PR values are systematically determined by the loading conditions of the bit in relation to the axis of material symmetry used to describe the rock formation.

In one specific implementation, the constitutive equations of linear elasticity are uniquely expressed through the application and use of MWD data to (i) populate the variables of the constitutive equations of linear elasticity and (ii) undertake an analysis of the constitutive equations to obtain measurements of near-wellbore mechanical rock properties of (a) Young's modulus of elasticity and (b) Poisson's ratio. Further, variations in the mechanical rock properties (e.g., YME and PR) are used to identify the nature and occurrence of mechanical boundaries or discontinuities in the subsurface such as fractures.

More specifically, a technique to determine near-wellbore mechanical rock properties, YME and PR, from MWD data may involve processing measurements of the weight-on-bit (WOB), torque-on-bit (TOB), Annular fluid pressure (Ap), angular bit speed (RPM) and components of motion describing the acceleration of the bit, including axial, and the rotary or tangential accelerations to (i) obtain sets of MWD data corresponding to known temporal and spatial positions along the borehole, (ii) calculate the forces acting on the rock formation in connection with the drilling apparatus and drilling fluids, (iii) calculate the displacements of the bit as it is accommodated by the deformation of the rock formation, (iv) inform the terms and loading conditions (variables) of a linear, elastic stress-strain relationship that describes the constitutive behavior of the rock formation in relation to the orientation of the well and (v) using the constitutive linear elastic equations as determined through the application and use of the MWD data to calculate the aforementioned mechanical rock properties, (YME and PR) and (vi) analyzing the YME and PR with respect to the axis of material symmetry in relation to the orientation of the well to identify the nature and occurrence of mechanical boundaries and discontinuities such as fractures and bedding planes among other things.

The present disclosure involves an innovative, new system, apparatus and method to specify, in general, the mechanical properties of a rock formation from an analysis of drilling vibrations generated by the cutting action of the bit and the deformation of the formation in response to the forces acting on the rock formation in connection with the drill bit and drilling fluid system while conducting drilling operations. Deformation may include elastic deformation, plastic deformation, and failure of the rock, which may be considered fracturing. Stated differently, aspects of the present disclosure involve obtaining information associated with the drilling of a borehole, while drilling, to identify mechanical rock properties of the formation being drilled. Such mechanical rock properties may be used, in some examples, to identify the presence of natural fractures or rock properties more or less susceptible to stimulation techniques. For example, knowing mechanical rock properties along a borehole or the presence of natural fractures along a borehole may be used to optimize hydraulic fracturing operations by focusing such fracturing on areas where it will be most effective, among other advantages. Mechanical rock properties may include elastic coefficients (e.g., the cij's), strength measurements such as initial yield strength, peak compressive strength, tensile strength YME, PR, shear modulus, bulk modulus, strain-hardening exponents, Thomsen coefficients and other mechanical rock properties.

The Nature of Elastic Coefficients as they Pertain to Fractures

The mechanics of drilling a well provide a natural, in-situ, means to measure the deformation of a rock formation and gather data suitable for determining mechanical rock properties, because the penetration of the drill bit is in and of itself accommodated by repeatedly fracturing the rock formation by using the bit to generate forces on the rock formation that are sufficient to overcome the failure strength of the rock, measurements of such in relation to the methods described here may be used in predictable ways to determine the presence of natural (in situ) fractures, fracture swarms (cluster of fractures), bedding planes, fault boundaries, and other information. In some instances, variations in mechanical rock properties are used to identify fractures, bedding planes and the like.

Elastic coefficients that describe a relatively large deformation in response to the forces acting on a rock formation indicate the rock formation is weaker and less rigid. Therefore mapping the spatial variations of the elastic coefficients provides information where there are zones of weakness in the rock formation. If the nature and occurrence of the zones of weakness in the rock formation as evidenced by changes in the elastic coefficients are localized in space or otherwise discrete relative to the surrounding elastic coefficients that would indicate the presence of a fracture or other mechanical discontinuity. Systematic changes in the spatial distribution of the elastic coefficients as are derived by the method are used to identify mechanical discontinuities and geological boundaries of rock formations such as bedding planes or faults that act to separate or offset rock formation with different rock properties, where the differences in rock properties are evidenced by the nature and distribution of the elastic coefficients.

As will be understood from the present disclosure, mechanical strength and deformation of the reservoir rock influences fracture creation, propagation and ability to maintain fracture permeability. FIGS. 1A-1B are simplified diagrams illustrating the difference between fractures induced in a relatively brittle rock formation, and which may include naturally occurring fractures, versus fractures induced in a relatively ductile rock formation, which may include fewer or no natural fractures, respectively. As illustrated in FIG. 1A, a horizontal section 10 of a borehole has been drilled through relatively brittle rock 12 and hydraulically fractured. In contrast, FIG. 1B illustrates a horizontal section 14 of a borehole drilled through relatively ductile rock 16 and hydraulically fractured. The fractures 18 created in the relatively brittle rock tend to penetrate deeper into the reservoir than the fracture 20 in ductile rock. Moreover, reservoir rock trending to the brittle end of the normal range tends to have higher initial production rates and lower decline rates. The techniques described in this disclosure will provide new information for selecting hydrocarbon bearing zones by differentiating between brittle rocks generally associated with larger fracture creation and better proppant support that is more permeable than ductile rock that produces smaller, less productive fractures which are prone to rapid compaction and closure and are less permeable. Similarly, techniques discussed herein may also identify areas where natural fractures may exist provide similar advantages as brittle rock. Generally speaking, as will be understood from the disclosure, various mechanical rock properties discussed herein will provide mechanisms whereby rock formations may be characterized, along the well bore, as to the relative brittleness or ductileness, or the relative susceptibility to stimulation techniques along the formation, which may include the identification of existing fractures or at least rock properties indicative of existing fractures

Data Acquisition Techniques

FIG. 2B is a diagram of a bottom hole assembly portion 20 of a drill string where the bottom hole assembly includes a drill bit 22, a mud motor 24, a bit sub 26 including various measurement components positioned between the drill bit and the mud motor, and sections of pipe 28 within a horizontal section 30 of a borehole 31, also referred to herein as a well bore. The vibration data used in the described methodologies may be recorded as close to the source (drill bit) as practical to avoid attenuation through the bottom hole assembly. One possible location for recording is directly behind the drill bit and ahead of the mud motor using the bit sub, although multiple bit subs may be used along the drill string for geophysical processing of the desired signal. Drilling a wellbore involves using a portion of the weight of the drill string, known as weight-on-bit, to push the drill bit into a formation 33. The rotating force on the drill bit, known as torque-on-bit, can come from the surface or from a mud motor close to the drill bit. When using a mud motor, drilling mud is pumped down the drill string until it encounters the power drive section of the mud motor where a portion of the mud pressure and flow is converted into a rotational force, which is mechanically coupled to the bit to thereby place rotational torque on the bit 22 to turn the bit. The rotational force on the bit can also be augmented by or come exclusively from mechanisms at the surface on the drilling rig.

The objective of the drilling process is to break the rock down into fragments that are small enough that they can be lifted and evacuated from the wellbore with drilling fluids in order to continue to accommodate the forwards motion of the bit. It should be noted that the action of the drill bit on a rock formation causes the fracturing of the rock formation along the borehole to drill the hole, and in the formation immediately adjacent the borehole. Moreover, the drill may encounter existing fractures 34 while drilling. Hydraulic fracturing, in contrast, is a process that occurs during the completion phase by injecting fluid into the borehole, typically with perforation clusters 22 in the casing, to initiate fractures 18/20 into the formation surrounding the bore hole, as illustrated in FIG. 1.

In the illustrated diagram, the bit sub 26 is shown between the bit and the mud motor. The bit sub is a cylindrical component that is operably coupled between the mud motor 24 and the drill bit 22 in a way that allows the mud motor to turn the bit. The bit sub provides housing, typically in a cylindrical shape, or mechanism to support various possible measurement components 36 including strain gauges, one or more accelerometers, pressure sensors, which may measure the pressure of the mud flow, temperature sensors which may measure the circulating temperature of the mud or other temperatures and which may be used to provide correction or offset of measurements or calculations that vary with temperature, gyroscopes which may be used to measure inclination and/or directional changes of the bit and string, and/or other components to measure or derive the information discussed herein.

In one example, as shown in FIGS. 2B and 2A, the strain gauges are mounted on the bit sub to determine torque on the bit and the weight on the bit (the force turning the bit and the force pushing the bit into the rock formation). Various possible ways of mounting the strain gauges, or combinations of strain gauges, are possible. Additionally, as shown in FIG. 2A, which is a representative front view of the bit 22, accelerometers are placed to measure axial, rotary, and/or lateral acceleration of the bit. Note, the bit axis is in the center of the circle, whereas axial acceleration may be measured somewhat offset from the axis depending on the placement of the accelerometer. Acceleration measurement may be accomplished by using one or more multi-axis accelerometers. The bit sub, or other such component, may also include a processor and memory to store computer executable instructions to implement various possible methodologies, and possibly preprocess data, as well as a power source which may be one or more batteries. Data storage, such as the memory or other data storage, is also provided to store the collected data. The measurement components, alone or in various possible combinations, may be provided in other locations of the drill string in the general proximity of the drill bit.

Geomechanics in Relation to the Cutting Action of the Bit

FIGS. 3A-3C are a sequence of diagrams illustrating a close up view of a cutter 32 portion of a drill bit in a borehole, slipping, sticking on a portion of rock, and then slipping loose when the forces on the bit are sufficient to overcome the rock causing the rock to fracture and the bit to rotate—collectively referred to as stick slip behavior. The rock deformation mode for a polycrystalline diamond compact bit (PDC) bit is shearing as opposed to a roller cone bit which is punching. Models that describe drilling behavior are in a large part informed on the mechanics of drilling with a roller cone bit and while these models have been extended for the application and use of PDC bits, they suffer uniquely from their inherent inability to reconcile the fundamentally different nature of rock deformation. As will be appreciated from this disclosure, the innovative, new techniques disclosed here seek to advance the application and use of PDC bits, as well as other bits, to characterize mechanical rock properties and in particular for the identification of the nature and occurrences of fractures. The figures describe the depth of cut in relationship to the area and displacement of the fractures created in response to the forces acting on the bit. The cutting action of a particular type of bit but should not be construed to limit the method in the use of other types of drill bits that generate acoustic emissions from rock failure in response to the forces acting on the geometry and configuration of the bit.

More specifically, as the bit turns, the interaction of the bit with the rock formation at any instant in time, produces a complex distribution of forces acting on the formation in connection with the bit and drilling fluid system (e.g., the mud motor) where the orientation and magnitudes of the forces acting on the rock formation are related to the configuration and geometry of the cutters on the bit. Generally speaking, drilling is not a smooth and consistent process. Instead, depending many things including the axial force on the bit, rotational torque on the bit, rock properties, and presence or absence or existing fractures, the bit cuts, gouges, spins, snags, and otherwise drills the borehole in a very complicated and varying fashion. In some instances, the complex distribution of forces acting on the formation is insufficient to initially overcome the strength of the rock formation in relation to the cutting action of the bit and the bit will stop rotating or stick.

As illustrated in FIGS. 3A and 3B, as the cutter begins to stick, the torque applied to the bit increases from a relatively steady value. As the forces, such as the illustrated torque-on-bit, applied to the bit change either through manual or automated interaction with the surface drilling apparatus or through the non-linear feedback of elastic energy stored within the drilling string or some combination of both, a new weight on the bit (WOB) and torque are delivered to the bit. These forces on one or more particular cutters will continue to load the rock elastically until such point (i) the rock begins to deform plastically and the deformation is concentrated along fracture planes and (ii) when those forces provides a sufficient distribution of forces to overcome the failure strength of the rock formation, the bit will turn and rock, often snapping loose, and drilling will continue. As shown in FIG. 3C, when the forces overcome the rock, the torque will dramatically drop to the relatively steady value, until the bit sticks again. Such a stick slip action may occur at varying frequencies and displacements and make happen one or more times per revolution of the bit per cutter on the bit, thereby resulting in many such cutting behaviors each revolution of the bit.

Regardless of whether and the extent of stick slip behavior is experienced, deformation and failure of the rock cause the bit to vibrate. During rock deformation and in particular when the bit overcomes the rock strength at failure, stored elastic strain energy is released in the form of acoustic emissions. In some instances, the bit is fracturing rock and may intersect fractures and existing mechanical discontinuities. In some instances, the bit may reactivate existing fractures, which may itself generate a distinct acoustical signal in the form of an induced bit vibration.

On the Processing of Drilling Vibrations

As introduced above, drill bits 22 typically include many cutters 32 arranged with a geometry and configuration designed to generate sufficient forces to overcome the failure strength of the rock formation 33 based on the nature of the rock formation expected to be encountered when drilling a well. During a single rotation of the bit, at least one, but typically many of the cutters will overcome the failure strength of the rock formation and produce a plethora of acoustical emissions related to the scraping, cutting, fracturing, and other interactions between the bit and the rock formation. The tool may include various possible mechanisms, including a reed switch or gyroscopes, that measure revolutions per minute and provide information of each rotation of the bit. There are as many as 30 cutter heads, so each rotation may cause hundreds of acoustic pulses.

Because of the stochastic nature of acoustic emissions or in relation to the nearly simultaneous initiation and propagation of multiple fractures at the cutting face, the implementation of the method discussed herein may use statistical methods and signal analysis tools. In one possible methodology, the Root Mean Squared (RMS) measurement technique determines the energy of the signal despite the shape of its waveform. This is important because the many simultaneous events will create complex waveforms with constructive and destructive interferences.

Given the stochastic nature of the acoustic emissions generated by the cutting action of the bit, it is expected that the more times the bit turns per unit time, the higher the rate of acoustic emissions. RMS levels obtained for a time window that underwent two revolutions of the bit would be expected to have higher RMS levels than those obtained during one revolution of the bit in the same time window, all other things being equal.

Measurements while drilling show that the bit speed can vary wildly and erratically during drilling operations. In some instances the bit can completely stick and then slip again as more force is gradually applied. If the stick-slip behavior of the drill bit is not accounted for in the geophysical signal processing, then the variations in the fracture measurements may be confused by with variations in the bit speed and not variations in the fracturing behavior of the rock. When the time window used to measure the RMS level of the signals that have been extracted from the drilling vibrations is normalized with respect to the bit speed, where the bit speed is recording using a gyro to sample the changes in position of the bit with respect to time or a magnet which is used to inform the position of the bit with respect to time, this normalized measurement is understood to provide a level of the acoustical emission activity generated by the fracturing of a rock formation in relation to the cutting action of the bit. The time window where data is gathered may be tied directly to bit rotation by having a time window set based on counts of the bit revolution. In other instances, the time window may be set, and the energy (from vibrations) may be normalized to account for some set turns of the bit, such as one turn of the bit.

The implementation of the technique as described uses signal processing techniques, such as Fourier transforms, bandpass filtering or other filtering, or combinations thereof, to calculate the motion of the bit and the forces on the bit from the amplitudes and frequencies of the acoustical signals recorded by the MWD apparatus (e.g. bit sub 26) that are generated in response to the cutting action of the bit (e.g., the drilling vibrations propagate up the drill string as acoustic waves sometimes referred to as collar waves or tool mode, where they are recorded as acoustical signals by the MWD apparatus). Forces on the bit may also be measured using strain gauges, in one possible implementation. The motion of the bit and the forces on the bit are used to populate a stress/strain relationship that allows the computation of mechanical rock properties. Mechanical rock properties may be analyzed relative to a baseline (such as an average over some wellbore distance) to identify locations where the mechanical properties of a rock formation change as the bit encounters a mechanical discontinuity or other geological discontinuity all other things being equal. In some instances, rock properties may be used to identify such locations through computations based on assumptions of the rock formation, and comparisons thereof, or otherwise.

Referring to FIG. 4, to account for the effects of the drilling efficiencies and the possibility of stick-slip behavior while conducting drilling operations, in one specific implementation, the method may use measurements of the drilling efficiencies, such as revolutions per minute (RPM) to normalize the geophysical signal processing of the drilling vibrations in order to compare the results of the signal processing along the length of the borehole in accordance with the methods provided. To get an RMS estimate, first square the amplitude of the data, here Axial acceleration. Then take the mean, either the length of time, T, or the number of samples in the frequency domain, N, where the number of samples may be band limited to between F1 and F2, but in general are understood to be from the first frequency in the Fast Fourier Transform (FFT) and the Nyquist frequency. Lastly, take the square root. It is understood through the normalization that the variations in the signal levels, such as the RMS levels, are now corrected in account of the changes in drilling efficiencies along the trajectory of the borehole.

As such, the RMS levels, obtained as shown above, for example, are related to changes in the mechanical rock properties while the drill intersects areas of differing properties. Stated differently, it may be helpful to compensate or normalize for changes in the rotational speed of the bit (RPM) by using a time window commensurate with the rotation of the bit. Such normalization may also be useful to correct the geophysical signal processing in accordance with drilling operations where the drill string is rotated from surface in conjunction with the mud motor turning the bit, as opposed to situations where the mud motor is operating but the string is sliding and not being rotated from the surface.

In a further elaboration of the method, the signal processing techniques systematically relate measurements of the forces acting on a rock formation in connection with the drill bit and drilling fluid system (stress) to the variations in the fracturing of a rock formation in response to the cutting action of the bit (strain) to obtain innovative, new stress-strain relationships where the application and use of the stress-strain relationships allow for the derivation of elastic coefficients for the stress-strain relationships. Relative variations in either one or a combination of the elastic coefficients may be used to identify the nature and occurrence of fractures, fracture swarms and other mechanical discontinuities and geological boundaries such as bedding planes and/or faults that offset or otherwise separate rock formations with different mechanical rock properties.

Accounting for the drilling efficiencies is also a consideration in the implementation of the approach, because the wear on the drill bit as the rock formation is drilled will change the configuration and geometry of the cutters on the drill bit. The mechanical wear of the drill bit will affect the distribution of the forces acting on the formation in connection with the drill bit and as the well is progressed, the manual or automated application of forces on the drilling apparatus change to account for the wear and tear of the bit.

The application and use of the stress-strain technique according the method employed herein is understood to normalize the effects of reduced drilling efficiencies caused by bit wear on the derivation of the mechanical rock properties because the forces used in the stress-strain relationship are obtained from the forces acting on the bit, where the forces needed to overcome the strength of the rock are increased in relation to the penetration or depth of cut of the worn bit tooth into the formation.

The acoustical emissions associated with the deformation and failure of the rock formation while drilling are generally too minute and/or too attenuated by the intervening rock to be detectable at the surface (which may be hundreds or thousands of feet above the borehole). Because of the amount of energy released is generally expected to be slight and of relatively high frequency, the radiated waves are best viewed when transmitted from the cutting face through the bit and bottom hole assembly where they propagate along the drill string through acoustically conductive steel as a direct tool arrival and contribute to the vibration of the drill string. The drilling vibrations can be recorded on instrumentation that is sensitive to their nature and presence. Stated differently, one aspect of the present disclosure involves a drilling tool assembly including sensors and processing electronics (e.g., the accelerometers and/or strain gauges in the bit sub 26 proximate the bit 22) that are positioned to detect and record the radiated waves from the drilling induced fracturing, which may further involve identification and/or characterization of existing mechanical discontinuities, such as fractures or geological boundaries, such as faults or bedding planes.

A Specific Data Logging Technique

In one specific implementation, a form of measurement while drilling (MWD) system or tool is employed. The MWD system uses sensors designed to measure vibration. The MWD system may also measure forces on the bit, and other parameters such as bit speed, which may be expressed as revolutions per minute, the fluid pressures, and temperature of the drilling mud or environment proximate the bit sub. The system may also include gyroscopes to obtain the orientation of the cutting face of the drill bit, in some implementations. In one specific embodiment, the MWD tool includes at least one receiver, accelerometers which may include strain gauges mounted on or proximate the bottom hole assembly to record the drilling vibrations and associated acoustic emissions. In some implementations, the MWD may further include electrical, mechanical, and/or other filtering mechanisms to processes the data to remove unwanted noise or to record the data without unwanted noise. In certain instances, stages of filtering may be applied both prior to recording, and after recording but prior to processing, to remove unwanted data, or as much as necessary or possible. In an alternative enablement, the signals may be transmitted to the surface for storage and processing. In some applications it may be desirable to process the acoustical signals, such as through the processor, on board the logging tool for transmission of the significantly data-reduced processed signal to surface in real time.

Once the drilling dynamics data is collected and processed, the results are correlated back to the measured depth of the well using precise measurements of the length of drill string components as they are lowered into the well. Gamma ray MWD and casing collar measurements can further be used to correlate the absolute location of data processed from the BHA collection point or points to determine a more reliable location of the bit in relation to the subsurface.

Noise Attenuation Techniques

Because of the amount of energy released is generally expected to be slight and of relatively high frequency, the radiated waves are best viewed when transmitted from the cutting face (with cutters 32) to the bit and bottom hole assembly, where they may further propagate and are known as the direct tool arrival or collar wave and contribute to the vibration of the drill string. The acoustic emissions may be measured by accelerometers, transducers, or other devices sensitive to particle motion.

Drilling induced vibrations that are generated by the interaction of the bit with the rock formation will have harmonic frequencies that are related to the rotational speed of the bit. Most of the harmonic vibrations are expected to be low frequency. The amplitudes of the acoustical emissions in the frequency ranges of the harmonic frequencies are usually much less than the amplitude of the harmonic drilling vibrations. In one implementation, the harmonic drilling vibrations are removed by a filter, such as a high pass filter or bandpass filter, which may be implemented in the bit sub processor, or may be applied to the stored data for later processing after download of the data from the bit sub memory, that passes signal frequencies that are higher than frequencies related to the harmonic drilling vibrations, and which may also eliminate frequencies above those possibly related to fracture characteristics. Other filter types and frequency characteristics will be possible depending on various factors including, but not limited to, the rotational speed of the bit, the type of bit, the rock characteristics, the positioning of sensors, mud motor characteristics, and other attributes.

Another consideration with respect to the variations in the amplitude and frequency of the acoustic emissions is the interference of the direct tool arrival by the generation and transmission of other wave modes that are also excited by the drilling operations including the wave modes excited by the release of energy from the fracturing of the rock formation.

In addition to the acoustical emissions created by the fracturing of a rock formation in response to the cutting action of the bit, the drilling vibrations that are being recorded will also generate wave modes that propagate though the formation surrounding the wellbore and the fluids within the wellbore that can interfere with the acoustical emissions that are related to the fracturing of a rock formation at the cutting face and bias the measurements of the mechanical rock properties. For some frequencies, the other wave modes will have higher amplitudes than the amplitudes of the acoustical emissions at that particular frequency. These other wave modes will result from the propagation of various guided waves in the drilling fluid between the bottom hole assembly and the wellbore such as Stonely waves, tube waves and direct fluid waves such as the fluid compressional waves. Other wave modes that can interfere with the direct tool arrival are surface waves that propagate and refract energy along the interface between the well bore and the fluid such as compressional head waves, surface compressional waves and shear body waves. These waves all have the potential to interfere with the propagation of energy from the direct tool arrival and if their presence is included in the signal processing could bias the calculations related to the amplitudes and frequencies used to describe the nature and occurrence of the fracturing.

To decrease the possibility of interference with the other wave modes, in one embodiment the sensor is mounted internally using a plug which will effectively isolate the sensor from the wave modes propagating through the formation and through the fluids. For example, in the case of accelerometers, the accelerometers are mounted within the bit sub (or other component). External waves that are carried by the formation and fluid in the annulus will therefore be mechanically isolated by the internal position of the sensor relative to the other waves. An internally mounted sensor will respond mainly to the vibrations related to direct tool arrival (the vibrations caused by the interaction of the bit with the formation).

Because the steel used in the construction of a bottom assembly has a significantly higher quality factor, Q of 10,000, that rock formations and fluids where Q ranges are typically from 1 to 100, the waves that propagate along the collar known as the collar wave or tool arrival will experience much less attenuation than signals recorded by the sensor that have travelled through the formation or the fluid media. Because of the low attenuation of the waves propagating along the steel drill collar relative to the waves propagating through the formations and drilling fluids, it could also be possible to naturally attenuate the various unwanted modes of propagation by placing the receiver at a distance that is far enough away from the bit to attenuate the other unwanted formation and fluid wave modes that will interfere with the direct tool arrival, but not so far as to lose the valuable high frequency information carried by the collar wave or tool wave that is needed to be recorded to calculate the size and displacement of the fracturing. The distances needed to attenuate the unwanted modes can be determined by expressions that relate the energy loss per cycle during transmission for a given quality factor. In this embodiment the location of the bottom hole assembly in relation to the location of the bit is placed at a distance behind the bit to achieve this attenuation. Alternatively or additionally, to avoid the interference of fluid and formation arrivals with the collar arrival the receiver may be coupled to the drill string at the surface of the well where the drill string has yet to enter the subsurface and borehole that contains the drilling fluids.

In another embodiment, a plurality of receivers are arranged along the string behind the drill bit to record the acoustical signals generated by the release of elastic energy at the cutting face. In one example, the receivers are spaced within centimeters or millimeters and may be placed in array mounted on a portion of the bit sub. The spacing of the receivers may form an array, with the spacing determined by the frequency range of the acoustical signals and the velocity of propagation in the steel or other material. Because the velocity of propagation in the steel is typically much faster than most formation and fluid velocities which control the nature of propagation of the interfering modes and is known with a high degree of certainty, spatial filters such as FK filters that pass signals that propagate at velocities that are consistent with wave transmission through steel and attenuate events with slower velocities. These filters can be used to separate the direct tool arrival waves from the other interfering waveforms. When using an array of receivers the signal-to-noise ratio can be further increased by using geophysical signal processing techniques to filter the data by stacking the signals over the array. The nature of the stacking depends on the configuration of the array of receivers and whether the receivers are deployed in a linear, bipolar or radial array.

In some instances the stacking can be used to isolate various modes that are propagated by the tool arrival, such as the compressional wave, transverse wave or quadrupole wave. These other tool modes can also be processed using geophysical signal processing techniques to determine the fracturing of a rock formation in relation to the cutting action of the bit.

The quadrupole is a direct tool mode that does not propagate direct tool arrivals above a cutoff frequency, where the cutoff frequency is related to the diameter and thickness of the steel. The amplitudes and frequencies of these other modes provides useful though otherwise band limited information that may be used in relation to the geophysical signal processing techniques in order to specify how the size and displacement of the fractures would be responsible for generating these other wave modes.

The analysis of the signals extracted from the drilling vibrations should not be limited to the case of the amplitudes and frequencies of the direct tool arrival. In another embodiment, the receiver array can be used to reject the direct tool arrival or collar wave and pass other arrivals related to other modes of transmission, where the modes of transmission may be through the rock formation or drilling fluid system based on their velocity of propagation and frequency content through the media. These other modes of propagation may be used in preference to the direct tool arrival when the other modes of wave propagation contain signals related to the fracturing of the formation in relation to the cutting action of the bit that are of interest to the variations of the RMS levels of the signals as described by the method.

The Application and Use of Microseismic Signal Processing Techniques in Relation to the Method

In a further elaboration of the method, drilling vibrations generated by the cutting action of the bit may be processed using geophysical signal processing techniques that are conventionally recognized as appropriate for the analysis of microearthquake source mechanisms. In one example, the depth of cut or penetration per revolution of the bit is obtained by using signal processing techniques to measure the sizes and displacements of fractures that result from the deformation and failure of the rock formation in relation to the cutting action of the bit, to estimate the zero-frequency level (ZFL) of the displacement spectra. The ZFL of the displacement is the static offset level here understood to represent the penetration of the bit. Under this consideration and in relation to the technique it provides high resolution motion of the bit that is taken to be the depth of cut or penetration of the bit for each revolution of the bit.

FIG. 5A illustrates a displacement spectra for one revolution of a bit. The graph of a displacement spectra has a y-axis of rotary displacement amplitude at a range of frequencies (x-axis) over which amplitude measurements are taken (based on the sampling frequency.) The displacement spectra is obtained by integrating the Fourier transform of the time domain measurement of the acceleration data (bit vibration data) twice in the frequency domain, or the angular velocity time series once where it is understood that the accelerations may be axial, lateral, rotary or a combination of those channels. Rotary displacement may be in the form of radians. The ZFL (zero frequency level) determined by the method is where the displacement spectra theoretically intersects the zero-frequency axis. The ZFL is directly proportional to the average displacements of the fractures initiated by the bit, and to the displacement of the bit per revolution. The high-frequency falloff in the rotary displacement spectra goes as ω⁻² which implies that the dislocations or rock failures caused by the cutting action of the drill bit can be represented by mathematical functions used to describe dislocations over a fault. This type of rock-failure behavior lends itself to a method to construct and analyze force-displacement or stress-strain diagrams in relation to rock deformation and rock failure using the ZFL to obtain the static dislocation and the accelerations as a force per unit mass.

Generally speaking, because high resolution displacement data is achievable, the system can detect relative changes of displacement per revolution of the bit or displacement per time or rate of penetration (ROP) and thereby determine when a fracture is encountered as the displacement will be greater relative to areas where fractures are not encountered. With respect to the method, the increase in displacement would indicate a change in the hardness of the formation, where hardness increases with displacement all other things being equal.

Observations of strong ground motion generated by earthquakes suggest that the time series recordings of the P- and S-wave signals can be reasonably and effectively treated as band limited white noise where the amplitudes and frequencies of the signals are controlled by the interaction of many smaller faults and fracture patches rupturing simultaneously and the band limitations in the absence of attenuation are related to the displacements and rupture dimensions of earthquake source. Therefore, observations involving the simultaneous occurrence of multiple acoustic emissions generated by fracturing at the face of the bit and then transmitted through steel where the signals are expected to undergo little attenuation indicates that the application and use of models that are generally used to describe the source mechanisms of an earthquake can be used to describe the aggregate sizes and displacements of many fractures being generated simultaneously from the repeated turning of the many cutters on the drill bit when overcoming the rock strength to accommodate the motion of the bit.

The estimation of the size and displacement of the fracturing generated by the cutting action or motion of the bit follows directly from the display of the data that is present in FIG. 5A. Here the model for representing the microearthquake source mechanism that estimate the size and displacement of the fracturing typically involves the application and use of a two parameter spectral model where here the two parameters used to describe the model are the zero-frequency level of the displacement spectra and the corner frequency. It should be appreciated that other source parameter models such as the RMS stress drop that also use the amplitude and frequencies of the signals to estimate the size and displacement of fracture and can also be considered and the use of a two parameter spectral model should not limit the scope of this method.

In the two parameter model of the microearthquake source mechanism, the displacement of the bit as determined by the displacements on the fractures is related to the zero frequency level of the displacement spectra. The ZFL of the displacement spectra is the same as the static offset. The low-frequency content may be modified by the receivers and electronics used to record the signals on the MWD assembly and the possible use of band pass filters in relation to the signal processing to eliminate the low frequency drilling harmonics. The technique has advantages because it (i) provides an estimate of the ZFL using band-limited data that may not always be reliable in the lower end of the amplitude spectra, (ii) provides a reliable technique to objectively select the ZFL without manual or visual biases, (iii) can be automated in a computer implemented fashion to handle the large volumes of data typically collected by an MWD apparatus as used to collect data in relation to the method.

Geophysical signal processing techniques that use relationships between the power spectral density of the displacement and velocity spectrum are used to calculate this value (ZFL) based on the functional mathematical representation of the earthquake source spectra as described by the two parameter spectral model of the earthquake source provide an objective technique to overcome the expected poor signal to noise ratios of the low frequency acoustical information.

Thus, ZFL is understood to represent the displacement or penetration of the bit. In one specific example, the ZFL is a measurement of the displacement (e.g., in millimeters or inches or angular displacements such as a radians) of the bit per turn of the bit. If the displacement is the axial displacement per revolution of the bit, then this displacement can be used to inform a depth of cut in terms of displacement per revolution. If time taken to make the one turn of the bit is used to describe the axial displacement or penetration of the bit, then it is taken that this may be used to determine rate of penetration. Typical depths of cut as estimated by the method range from 0.01 inches per revolution to 0.1 inches per revolution

The size or the radius of the fractures can also be determined by the frequency content of the acoustical signal. Signals with higher frequency content generally correspond to smaller fracture areas. This characteristic of the frequency spectra that is used to determine the size of the event is typically referred to as the corner frequency (FIG. 5A). There is a linear, relationship between the corner frequency of the displacement spectra and the size of the event. Because of attenuation of the signal during the transmission, estimates of the corner frequency may be compromised if the sensor is placed too far from the bit. Thus, having the sensors as close to the bit as possible, or at least not so far that attenuation is significant, may be a consideration for some implementations discussed herein. Limitations in the sensor and the recording electronics and the filters employed to extract the signals may also limit the useable bandwidth.

Geophysical signal processing techniques that use relationships between the power spectral density of the displacement and power spectral density of the velocity based on the functional mathematical representation of the earthquake source spectra as is described by the two parameter model of the earthquake source that can be used to calculate the corner frequency. A relatively high corner frequency may represent unfractured rock whereas relatively lower corner frequency may indicate the presence of a fracture.

When the bandwidth is limited, the minimum fracture size detected will be set to a threshold based on a cutoff frequency. Any fracturing of a rock formation below this threshold will not contribute to the determination of the fracturing of the rock formation and the deformation of the rock formation. Thus, in the application of the method, some rock formations may undergo fracturing in response to the cutting action of the bit where the fracturing is not detected by the method employed here, either because the frequencies of which to fracture energy occur are too high or the signal to noise ratio is too low. The occurrence of these scenarios would suggest that the bit is not effectively penetrating the formation and that the rock properties derived in this instance would represent the limiting rock properties and limit the specification of the deformation of the rock in response to forces generated by the emplacement of hydraulic fractures where any such deformation in relation to the cutting action of the bit would be taken as insufficient to do so.

In order to make meaningful comparisons of the fracture sizes and displacements along the well bore, the time window that is used to process the signals extracted from the drilling vibrations is based on the bit speed, in one specific implementation. One way to normalize the time window relative to the bit speed is to specify the time window according to the time needed to make one revolution of the bit. If the bit speed were 120 RPM then the time window would be 500 ms, while if the bit speed were 60 RPM then the time window would need to be 1000 ms to obtain an equivalent measurement of the fracturing generated in response to the cutting action of the bit. Another way to normalize the time window to account for variations in the bit speed along the wellbore would be to normalize the temporal frequency of the spectral density by converting from cycles per second to cycles per revolution. At a bit speed of 120 RPM, the bit would make two revolutions in a second and therefore at cycle per second would be normalized to two cycles per revolution, while a bit speed of 60 RPM would make one revolution in a second therefore one cycle per second would be normalized to one cycle per revolution. In one embodiment, the BHA is instrumented to measure the bit speed at sufficient resolution to specify either a time window equivalent to one rotation of the bit or normalize the spectral density by converting the temporal frequency to cycles per revolution.

As such, the spatial variations in the measurements as obtained through a combination of one or more of the measurements such as the RMS acceleration or the sizes and displacements of the fractures are understood to correspond to the nature and occurrence of deformation and failure in relation to the cutting action of the bit and as such are taken to represent the spatial variations in mechanical rock properties. So, an increasing value of ZFL, relative to a baseline, represents the intersection of the bit with zone of mechanically weaker rock that may be determined to be a fracture, swarm of fractures (e.g., fracture 34) in accordance with for example a stress-strain relationship. Therefore, the geophysical signal processing techniques employed by the method may involve statistical descriptions such as measuring the RMS level of the acoustical emissions or the application and use of geophysical signal processing techniques that are generally recognized as appropriate for the analysis of microearthquake source mechanisms to describe the deformation and failure of a rock formation in relation to the cutting action of the bit, where the spatial variations in the measurements (e.g., changes in ZFL or corner frequencies) relative to some average or baseline level are used to identify mechanical discontinuities or geological formations as they are encountered or crossed by the bit.

In another embodiment, referring now to FIG. 6, the geophysical signal processing methods may obtain measurements related to the instantaneous change in the inclination of the bit relative to the average direction of the bit to describe a mechanical discontinuity or geological boundary. Instantaneous change in inclination of the drill bit relative to a long-term average inclination can be obtained using a sensor or an array of sensors configured to record and extract acoustical signals in relation to the independent spatial axes of the drilling vibrations and will show the orientation of natural fractures drilled obliquely. The intensity of the inclination will be proportional to the angle of intersection. The magnitude of the deflection of the inclination of the drill bit relative to a long-term trend in the direction of the drill bit depends on the orientation of the mechanical discontinuity or geological boundary with respect to the cutting face of the drill bit and therefore the deflection is understood to indicate a change in the mechanical rock properties. So, for example, with a multi-axis accelerometer measuring axial acceleration and lateral or rotary acceleration of the bit, the ratio of the axial and lateral or rotary displacements may be treated as an inclination, as illustrated in the respective inclination logs 60 (prior to the fracture, while the bit 22 intersects the existing fracture 34). The lateral or rotary acceleration will increase relatively when the bit is deflected from axial movement, such as when the bit encounters an angled discontinuity transverse the borehole.

The Concept of Specific Energy in Relation to the Method

In addition to acoustical information that informs the displacement of the bit, aspects of the present disclosure may further involve force information. While conducting drilling operations the energy defined as the energy needed to remove a volume of rock is useful to describe the efficiency of the drilling operation. The specific energy is a term that describes the minimum amount of work needed to remove a certain volume of rock. Descriptions of the specific energy are useful to understand the variations in rate of penetration or the depth of cut in relation to the forces acting on the bit. Unlike conventional methods, the specific energy here may be calculated from the aforementioned acoustical processing techniques in and of themselves or in conjunction with measurements of torque-on-bit and/or force on bit.

The size and displacement of the fractures that form in response to the cutting action of the bit controls depth of cut into a rock formation and subsequently the rate of penetration of the drill bit through a rock formation. In a method described here, the work done per the volume of rock removed is determined by the displacement of the fractures as multiplied by the area of the bit. The work is computed from the forces acting on the bit multiplied by the displacement of the bit. As discussed above, the displacements of the bit are obtained through the analysis of the drilling vibration using geophysical signal processing techniques that are appropriate for the analysis of microearthquake source mechanisms as is provided by the above discussed method.

There are two components of the specific energy: one that is normal to the bit and another that is tangential to the bit or rotary energy as shown in FIG. 5B. Angular acceleration is the rate of change of angular velocity. The rotary displacement is obtained by dividing the angular velocity spectra by 2πf or by dividing the angular or rotary acceleration by (2πf)². Because the rotary specific energy is proportionally related to the torque per unit of displacement, it provides a measure of the rock strength or the minimum energy needed to drill and the application and use of the rotary specific energy in this regard takes the form of a stress-strain relationship (FIG. 8A and others).

In the method provided here, the displacement of the bit as it relates to the fracturing of the formation is evidenced through the RMS level of the measurements or the displacement of the bit is evidenced by the displacement on the fractures as provided through the microearthquake source parameter measurements and the volume of rock removed is proportional to the fracture areas and the fracture displacements are averaged over all of the fractures provided in the period of time processed for the time period that is analyzed for example in a single turn of the bit. In one embodiment of the method, the volume of rock excavated is specified by the product of the aggregate area of the fractures and the average displacement of the fractures removed by one turn of the bit.

For a bit that is turning at a rate of 120 RPM this would involve using a time window of 500 ms to the average area and average the displacement for one revolution of the drill bit. It should be appreciated that in instances when the rock properties are varying slowly or the MWD measurements are updated at rates less than the time period of one bit revolution the method is not limited to periods that are specified by the turning rate of the bit.

The rotational specific energy is the torque over the average fracture displacement per revolution as obtained by the analysis of the vibration data disclosed by this technique, provides a novel, innovative stress-strain relationship (FIG. 8A and others). For a given rock type, this is expected to be a linear relationship where the slope of the line is related to an elastic coefficient describing the strength of the rock (e.g., specified rock strength 1 (70), specified rock strength 2 (72), and specified rock strength 3 (74)). Because the measurements are obtained while drilling, the strength of the rock actually determined by and is intrinsically related to the fracturing in response to the cutting action of the bit. When using a PDC bit, for example, the strength of the rock determined from this stress strain relationship would be directly related to the shear strength of the rock. Stress-strain measurements obtained using the method disclosed can be used to characterize the elastic coefficients of the rock formation.

Innovative Stress Strain Relationships as are Provided by the Method

The stress strain relationships employed by the method are populated from measurements taken while drilling. In one instance, the strain is understood to be related to the depth of cut or the penetration per revolution of the bit which is determined by differencing the spatial location at two instances in time versus the number of revolutions taken for the bit to travel that distance. In accordance with the methods set out here in, where (i) the drilling vibrations are understood to represent the deformation and failure of a rock formation in response to cutting action of the bit in order to accommodate the forward motion of the bit through a rock formation, and (ii) the drilling vibrations as processed through the signal processing techniques to evidence the motion of the bit per turn of the bit are understood to represent strain and (iii) the RMS acceleration as obtained through the geophysical signal processing techniques to evidence the forces acting on the bit per turn of the bit inare used to populate the variable in relation to a stress-strain constitutive equation. In a further elaboration, these strain measurements can be related to the forces such as the weight-on-bit and torque-on-bit acting on the formation in connection with the bit and drilling fluid system to provide a diagram of a general stress-strain relationship. By relating the orientation and magnitude of the stress with respect to the orientation and magnitudes of the strain, where the orientations and magnitudes of the stress and strain are related to the geographical coordinates of the well, multiple stress-strain relationships can be established to determine the elastic coefficients of a rock formation. Where the well is drilled perpendicular to the maximum horizontal compressive stress, it is understood that these stress-strain relationships are expressed in the principle axes.

In general, the deformation of a homogeneous isotropic rock formation can be specified by two elastic coefficients. As the complexity of the rock formation increases through the presence of mechanical discontinuities and geological boundaries, the number of elastic coefficients needed to fully describe the deformation of a rock formation in response to the forces acting on the formation in general increases.

The simplest elastic coefficient would be to relate the WOB to the strain generated by the cutting action of the bit where the strain is the displacement in relation to a length made by one turn of the bit where the time windows used for the geophysical signal processing techniques are related to the bit speed. In one preferred embodiment of the method, the WOB would be obtained from the RMS acceleration, where the component of acceleration is oriented parallel to the borehole and as is illustrated in FIG. 10B.

For a transverse isotropic (TI) media, the stress strain relationship constitutive equation is illustrated in FIG. 8A (using MWD data) with attendant PR intercept and YME slope information illustrated in the curves of FIG. 8B.

Where:

E is Young's Modulus of Elasticity (YME),

ν is Poisson's ratio (PR)

σ1=TOB

σ3=WOB

ε3=Axial displacement spectra ZFL

This stress-strain relationship would in general be proportional to a Young's modulus of elasticity where YME is determined parallel to the direction of drilling. The method allows for the determination of elastic coefficients for transverse isotropic elastic media, where the assumption of transverse isotropic elasticity is understood to be reasonable approximation to describe the rock deformation by obtaining the stress-strain relationship to the orientation geological boundaries with respect to the inclination of the bit of the well being drilled. This could be accomplished using FIG. 8A or FIG. 10C where the elastic coefficients as described are proportional to Young's Modulus and the Poisson's ratio of a rock formation.

In another implementation more particularly shown in FIGS. 9A and 9B, the rock strength is specified by the stress-strain relationships based on the orientation of the well in relation to the bedding planes of the hydrocarbon bearing formation and the orientation of the well with respect to the principle axes of tectonic stress or the state-of-stress acting on the rock formation. This method allows for the determination of elastic coefficients for a transverse isotropic (TI) elastic media, where the assumption of transverse elasticity is taken as a reasonable approximation to describe the rock deformation by obtaining the stress-strain relationship to the orientation geological boundaries with respect to the inclination of the bit of the well being drilled. This could be accomplished using the above description where the elastic coefficients as described are proportional to Young's Modulus and the Poisson's ratio of a rock formation.

The loading conditions on a rock formation are given by the forces acting on the rock formation in connection with the drilling apparatus (e.g., weigh on bit and/or torque-on-bit) and drilling fluid system (e.g., annular pressure) and the deformation of the rock formation are described by the displacements of the bit (e.g., axial and lateral or rotary displacements as measured by accelerometers in the bit sub). The constitutive equations are described with respect to TI media. As shown in FIG. 9A-9B, TI media involves a layering of media (sometimes referred to as “layercake”). Transverse isotropic media involves a layering of media that is normal to a plane of isotropy—meaning the media is relatively uniform around the axis of symmetry. FIG. 9A illustrates a case where the axis 90 of material symmetry is parallel to the borehole 94 (and parallel to the bit axis 96 drilling the borehole 94). FIG. 9B illustrates a case where the axis of material symmetry 98 is perpendicular to the borehole 100 (and parallel to the bit axis 102 drilling the borehole 100). In one possible implementation, it is assumed that the media is either vertically transverse isotropic (VTI) or it is horizontally transverse isotropic (HTI). VTI is a case where the axis of symmetry is vertically oriented (layers are horizontal with respect to the free surface). HTI is a case where the axis of symmetry is horizontally oriented (where the anisotropy is understood to involve a layering of the formation that is vertical with respect to the free surface and further where it is understood that fractures represent a case of vertical layering with respect to the free surface). Thus, in the case of a vertical well, FIG. 9A illustrates vertical transverse isotropy (the vertical borehole is parallel to a vertical axis of isotropic layer symmetry) and also illustrates, in the case of a horizontal well, horizontal transverse isotropy (the horizontal borehole is parallel to a horizontal axis of isotropic layer symmetry). In contrast, in the case of a horizontal well, FIG. 9B illustrates vertical transverse isotropy (the horizontal borehole is perpendicular a vertical axis of isotropic layer symmetry) and also illustrates, in the case of a vertical well, horizontal transverse isotropy (the vertical borehole is parallel a horizontal axis of isotropic layer symmetry).

FIGS. 9A and 9B also illustrate the variables for the constitutive stress-strain equation (FIG. 8A and FIG. 10C) of a TI media. The stress strain variables are populated based on either (i) force measurements (e.g., from strain gauges) and/or (ii) the acceleration measurements (e.g., from accelerometers) or populated only from acceleration data. More specifically, in the first case, the stress-strain variables are populated with data related to WOB, TOB or axial or lateral or rotary displacements. In the second case, the stress strain variables the forces are populated with axial and lateral or rotary acceleration data, and the strain from the axial or lateral or rotary displacements (which may be available by integrating the axial or lateral or rotary acceleration spectra twice in the frequency domain). Generally speaking, Poisson's ratio (and/or Young's modulus) is computed from stress-strain constitutive equations populated from measurements taken while drilling along a borehole under the assumption that the transverse axis of material symmetry is parallel to the borehole (FIG. 9A) and under the assumption that the transverse axis of material symmetry is perpendicular to the borehole (FIG. 9B). FIG. 8B is a graphical representation of these constitutive stress-strain equations and illustrates a curve (solid line) where the axis of the TI material symmetry is parallel to the well axis (FIG. 9A) and a curve (dashed line) where the axis of the TI material symmetry is perpendicular to the well axis (FIG. 9B).

When the rock formation (media) is isotropic, the two equations will generate values of PR or YME that generally track each other, meaning that under either assumptive case, the drill bit will react similarly as it drills through relatively uniform rock and therefore the two calculations of PR and YME, while different, will nonetheless track each other. If, however, the borehole intersects discontinuities and or more the media is anisotropic, the computations of PR and/or YME will no longer track each other. The stress-strain relationships may be recast into other equivalent forms, and the particular arrangement shown in FIGS. 8A, 9A-9B, and 10C are taken for the sake of conveniences and should not be considered limiting. In a particularly useful manner as shown in FIG. 11, the constitutive equations could be re-arranged so that the Young's modulus term is determined using the slope of a linear relationship if the equations were cast in terms of the ratio of the forces acting parallel to the axis of symmetry to the forces acting perpendicular to the axis of material symmetry.

In certain instances, such as drilling under high confining annular pressure (AP), the confining pressure may be more accurately described using the AP instead of WOB. In such a case, AP would be used in place of WOB for the cases illustrated in FIGS. 9A and 9B, and possibly others.

Some of the stress strain relationships are set out in terms of bi-axial loading accounting for weight-on-bit and torque-on-bit. However, it should be understood that an implementation accounting for tri-axial loading may nonetheless use and account for WOB and TOB or the RMS acceleration.

Thus, in the application of the stress-strain relationships for a layered media, where the bedding planes are horizontal a vertical well and a horizontal well are drilled through the same rock formation where:

-   -   1. Separate stress-strain relationships that describe the         elastic coefficients when the bit is cutting perpendicular to         the geological boundaries such as the bedding planes and when         the bit is cutting parallel to the geological boundaries would         enable at least four elastic coefficients to be determined which         can be used to describe in general, stress-strain relationships         of transverse isotropic rock formations.     -   2. Separate stress-strain relationships that describe the         elastic coefficients when the bit is cutting parallel to the         direction of maximum horizontal compressive stress and parallel         to the maximum vertical compressive stress would enable at least         four elastic coefficients to be determined.

Referring to FIG. 10C, under the two assumptions (TI axis parallel to borehole and TI axis perpendicular to borehole (FIGS. 9A, 10A and 9B, 10B)), there are two values for PR and two values for YME, with PR and YME being elastic coefficients.

In practical application, the rock formation may not be actually horizontal or vertical, but may be tilted. When drilling through titled media, particularly media with small deviated angles of less than 30 degrees, given the expression of the trigonometric variables in relation to a rotation of the principle axis, the equations still produce useful results for PR and YME and may be used to inform the variation in the mechanical rock properties and in particular inform the location of fractures where the variation of the mechanical rock properties is predicted to do so as is described below.

On the Calculation of the Elastic Coefficients YME and PR Through a Mechanical Rock Property Analysis (MRPA) of the Method

FIG. 12 is a diagram illustrating linear stress strain relationships with curves fit to data pairs corresponding to different locations along a bore hole. The slopes of the fit lines in these examples relate to YME in different locations along a bore hole. To obtain data for populating the constitutive equation (or equations), various geophysical data processing techniques are involved, where:

-   -   1. Sampling the MWD measurements at a sufficiently high         frequency to resolve a small degree of mechanical variability         (which may correspond to the nature and occurrence of a discrete         fracture a couple of mm wide), the measurements of the forces         acting on a rock formation in connection with a drilling bit,         (TOB and WOB) and fluid system pressure or AP, the angular speed         of the bit expressed in revolutions per unit of time (RPM) and         the 3-components of motion that represent the acceleration of         the bit and where each of the measurements taken in time         corresponds to a discrete position of the bit along the length         wellbore (the MWD data)     -   2. Processing accelerations may be processed using geophysical         signal processing techniques to obtain (i) the average lateral         or rotary and vertical displacement of the bit that correspond         to a single revolution or turn of the bit and (ii) Root Mean         Squared (RMS) amplitude of the acceleration     -   3. Calculating the RMS averages of the TOB and WOB over time         windows that correspond to a single bit turn. Typical         penetration rates are usually 0.02 in/rev, while drilling at 240         rpm, would result in a single turn of the bit every 250 ms which         if sampled at 1 kHz would provide sufficient data to be able to         identify when the bit encounters a single-discrete fracture     -   4. Obtaining data pairs, which collectively define the linear         stress strain relationship, from MWD data using relationships         between the RMS acceleration and displacement of the bit that         are appropriate for the loading conditions and motion of the bit         in relation to the axis of symmetry for the constitutive         equations used to describe the rock formation where it is         understood that the loading conditions are determined with         respect to the orientation of the drilling well.     -   5. Using curve fitting techniques to estimate the two parameters         needed to describe a line, the Slope and Intercept and         statistical descriptions of the variations of the two elastic         parameters YME and PR with respect to each of the linear         clusters as were identified with respect to the band-limited MWD         data that was used to generate the data pairs     -   6. Identifying where the distributions of the MWD data pairs         form spatially or temporally consistent clusters or otherwise a         locus of adjacent points that (i) can be described through the         application and use of curve fitting techniques in terms of         relationships that are linear and (ii) such that each of the         linear relationships can be used to determine the parameters of         line such as the slope and intercept in relation to the         mechanical rock properties YME and PR.

In one possible example, a set of data pairs are generated for the equation of FIG. 9A (10A) and/or the FIG. 9B (10B). The data pairs are generated along the length of a well bore. In one example, a data pair comprises (y, x) and the set of data pairs may be used to generate the linear stress strain relationships illustrated in graphical form in FIG. 12 for each of the cases shown in FIGS. 9A and 9B.

In one embodiment, the variables of the equations are populated using (i) the RMS acceleration data, which is used to describe the forces acting on the formation in connection with the bit, and (ii) the ZFL of the displacement spectrum is used to describe the motion on the bit where the motion understood to be the strain experienced by the rock formation and where the orientations of strain are described by the bit displacements according to the cutting direction of the bit and the orientation of the borehole in relation to the orientation of the axis of material symmetry as is shown in FIGS. 10A-10B. Thus, the acoustical acceleration measurements may be use to specify the forces acting on the formation. The approach essentially conforms with Newton's second law and balance of forces.

Special Cases to Consider

As shown, the MWD parameters may be expressed as a function of the frequency content. In particular, where the frequency content is limited through the application and use of a band-pass filter to generate band-limited MWD data pairs to populate the terms and conditions of the constitutive equations of linear elasticity, using band-limited MWD data to form frequency-dependent, data pairs corresponding to data that generally describes the stresses acting on the material and the deformation of the material to populate the terms and conditions of the constitutive equations of linear elasticity. Using geophysical signal processing techniques to obtain zero-frequency levels (ZFL) of the displacement spectra where the ZFL corresponds to depth of cut per revolution of the bit or the depth of cut per unit of time or the rate of penetration (ROP) in relation to the width of the bandpass filter used to window the MWD data (the band-limited MWD data). That is the ZFL is determined from a specified range of frequencies or is otherwise calculated from bandlimited data.

The two parameters YME and PR and the statistical descriptions of the variations of the two parameters as can be obtained through the curve fitting technique for each liner cluster can be used to diagnose the drilling conditions for each of the depths as a function of the frequency used to form the data pairs.

Fracture Identification from the Elastic Coefficients

It is usually not known a priori what the appropriate angular relationships between the axis of symmetry of the constitutive elastic equations used to describe the rock formation and the axis of symmetry of the wellbore are. In practice, most horizontal wells are drilled parallel to bedding and most vertical wells are drilled perpendicular to bedding. Further, in basins, most natural fractures are vertical, and thus most horizontal wells are drilled perpendicular to fractures and most vertical wells drill parallel to fractures.

Referring again to FIGS. 9A and 9B, the MRPA technique is used to populate the terms and conditions of the two assumptive cases where the axis of symmetry of the constitutive media used to describe the rock formation are examined for the cases when (i) the axis of material symmetry is parallel to the axis of drilling and (FIG. 9A) (ii) the axis of material symmetry is perpendicular to the axis of drilling (FIG. 9B). In an isotropic media, the determination of YME and PR from the MRPA analysis for the two constitutive stress strain relationships will result in equally or closely spaced values of YME and PR that tend to track each other in space and time. In anisotropic media, the two computations of YME and PR as provided by the MRPA will deviate in predictable ways that can be used to identify the nature and occurrence fractures in relation to (i) the differences in the elastic coefficients as specified by the type of anisotropy, either HTI or VTI, that is encountered relative to the orientation of the drilling well and (ii) the reduction in strength of the rock as is provided by variations in HTI YME and the VTI YME.

In the case of a horizontal well when the constitutive equations describe a variation in the end members situations (i) the rock may be understood to be fractured when PR VTI is lower than PR HTI and the magnitude of the fracturing is related to the decrease in the calculation of YME. Typically, in practice, fractured zones can be discerned in horizontal wells when the media VTI PR is lower than the HTI PR. The difference in the values of the elastic coefficients YME and PR between the HTI and VTI represent two cases of a cross-over curve or an anisotropic cross-over curve.

In the case of a vertical well, the rock is understood to be fractured when HTI PR is less than VTI PR (e.g., FIG. 13, discussed in more detail below). An objective way to determine these cross-over relationships (PR or YME for the two assumptive cases) is by calculating a series in time or a series in depth (the logs) of the YME or the PR for both the VTI and HTI solutions. Time or depth may be correlated to length along the well bore from which the measurements were made.

Mechanical rock property logs (the “Logs”) as calculated using the equations set out by the method can be smoothed by averaging the PR and YME values by (i) using the statistics in relation to the goodness of fit to the curve where one such statistic is known as the R statistic as provided by standard least-squares linear curve fitters such as LINEST in Excel, to filter out data with poor statistical evidence for a linear relationship between the data pairs or (ii) weigh the values of the curve at a particular location. Using a curve smoothing technique that would average the data over a time window where the length of the time window corresponded to the variation of the data and then resampling the time window in either time or depth according to its position in the subsurface. Processing logs using smoothing techniques can improve the ability to identify the relative variations in the elastic coefficients.

Further processing the Logs by subtracting the mean value from the HTI and VTI YME Logs where the mean value is a running average of the data along the Log where the length of the running average is related to:

-   -   (1) the measure of the smallest length of mechanical anisotropy         variation of interest which in terms of the practice and         application of the method may involve spatial distances as small         as 1 inch for the sample rate and frequency content that is         afforded by the application of use of state-of-the-art MWD in         relation to the method         -   AND     -   (2) of sufficient resolution to document the nature and         occurrence of changes in rock properties needed to identify         fractures at the resolution need for the commercial exploitation         of commercial hydrocarbons from unconventional reservoirs.

The subtraction of the mean value from the HTI and VTI PR Logs provides a baseline from which to compare the variations between the curves in ways that can be used to identify the locations of fractures. In practice the mean values of the Logs can be calculated as the average of the data values in the Logs over a certain time or distance specified by the spatial position in the subsurface along the length of the wellbore from which the measurements were made. When this mean value is subtracted from the Logs it provides the mean-subtracted Logs from which to make it convenient to obtain relative comparisons (e.g., PR HTI to PR VTI and/or YME VTI to YME HTI).

Taking differences of the mean-subtracted YME HTI Log and the mean-subtracted YME VTI Log and taking the differences of the mean-subtracted PR HTI log and the mean-subtracted PR VTI log provides an identification of the type of rock anisotropy based on the orientation of the well in relation to the axis of material symmetry. The differences in these Logs as evidenced by the behavior of the elastic coefficients in relation to the orientation of media symmetry with respect to the orientation of the wellbore can be can be understood in predictable ways to describe the location of a zone of weakness in relation to the method understood to be a fracture. For the cases presented here, these relationships provide a predictable way to identify fractured rock formations, among other advantages

Specific Case 1: Vertical Well, Horizontal TI Media

In typical VTI media, the vertical PR (FIG. 9A—TI axis parallel to borehole) typically has lower values than horizontal PR (FIG. 9B—TI axis perpendicular to borehole). That is, the material is more compliant to a load that is applied perpendicular to the axis of material symmetry than to a load applied parallel to the axis of material symmetry. Loading the rock formation in the same direction as the axis of material symmetry will result in less horizontal deformation, decreased horizontal compliance and/or higher ratios of horizontal to vertical stiffness. Conversely, loading the rock formation perpendicular to the axis of material symmetry, when the material symmetry is governed by fractures, will have higher compliances, higher PR and lower YME.

Stated differently, in a vertical well when the media is layercake (the axis of media symmetry is parallel the borehole), the weight-on-bit and axial displacements are typically parallel to the axis of material symmetry. In this case the VTI PR is typically less than HTI PR. Therefore, detection of zones where PR VTI is greater than PR HTI implies that the material behavior under the loading conditions of the bit is stiffer or less compliant in the horizontal direction as opposed to the vertical direction. Referring to FIG. 13, when a vertical well encounters HTI media (in a formation expected to be VTI), the mean-subtracted PR HTI log becomes less than the mean-subtracted PR VTI log, and a Fracture ID flag may be generated. This flag means the ratio of the displacement parallel to the axis of material symmetry increases relative to the displacement perpendicular to the axis of material symmetry. This is evidenced through the application and use of MWD data to describe the constitutive behavior of a rock formation and demonstrate an increase in PR VTI relative to PR HTI. This increase in PR evidences the presence of vertical fractures or a vertically fractured rock formation when drilling a vertical well. This logic is generally true in most circumstances, because most vertical wells are drilled perpendicular to the bedding planes of the rock formation where the well bore axis is parallel to the axis of material symmetry and so the presence of high VTI PR values is probably not related to vertical bedding planes. This behavior of a vertically drilling well encountering a set of vertical fractures can be further corroborated by a similar examination of the differences between the VTI YME and the HTI YME. Because the torque-on-bit is acting as a body force parallel to the axis of material symmetry as would be expected in the case of a vertically drilling well encountering a set of vertical fractures where the fractures control the axis of material symmetry, the VTI YME decreases and the HTI YME increases.

When the YME VTI log crosses or decreases relative the YME HTI log for a horizontal well then it is likely a zone of weakness that is more compliant in the direction of loading that is parallel to the WOB or perpendicular to the axis of material symmetry has been detected which would again be consistent with a zone of vertical fractures. So, in some instances, PR crossover would generate a flag, YME crossover would generate a flag, and in some instances the presence of both flags indicate a fracture. Moreover, in practice, a threshold may be applied that would need to be met before generating a flag. In one example, for PR data, a distribution curve may be generated for all positive crossovers, and only crossovers exceeding 68th or 90th or 95th percentile may be flagged. Other thresholds or data technique may also be used to eliminate data points that may be attributable to noise.

When the relationship between the HTI and VTI elastic constants returns to a pre-crossover relationship for a vertical well in an HTI media, the crossover Fracture ID flag is set back to zero. The flag may remain set, however, for as long as the log indicates. So, as shown in the example of FIG. 13, there are six sections of the 100 foot illustrated borehole where fractures are identified over several feet for each section.

Specific Case 2: Horizontal Well, Horizontal TI Media

FIG. 14A illustrates the two YME computations (FIG. 10A and FIG. 10B) over about 100 feet of horizontal well. FIG. 14B illustrates PR ratio computations (also FIG. 10A and FIG. 10B) over the same 100 fee of horizontal well. The YME and PR computations are based on data obtained while drilling. In this case, only the geophysical processing of the acoustical data as obtained from the accelerometers are used to populate the equations, so direct measurements of WOB or TOB are not provide in this example, but they could be included in another application and use of the equations as provided by the method. In the example illustrated in FIGS. 14A and 14B, a fracture ID 140 may be generated where shown. In this case, effectively the opposite behavior than was described for a vertical well drilling parallel to the axis of symmetry or VTI or layercake, is expected because now, in the case of a horizontal well in TI media, the axis of drilling is perpendicular to the axis of media symmetry. In this case, the loading conditions (torque-on-bit) are parallel to the axis of material symmetry conditions and rotational displacements (revolutions of the bit) are also parallel to the axis of material symmetry and HTI calculation will result in a lower HTI PR.

In the case of a horizontal well encountering a vertical fracture, the relationship is expected to be similar to a vertical well drilling in a VTI media—the PR HTI will decrease relative to the mean-subtracted PR VTI. More specifically, in the event with the mean-subtracted VTI PR for a horizontal well drilling in a layercake media is found to be less than the mean-subtracted HTI PR, then it is likely that the axis of material symmetry relative to the orientation of the drilling well and loading conditions of the bit has been rotated by 90 degrees. This can occur when the axis material symmetry is defined by a set of vertical fractures being intersected by a horizontal well. Thus, as shown, for example in FIGS. 14A and 14B, a fracture flag may be identified in the identified areas, as well as possibly other areas.

Likewise, the mean subtracted VTI YME will increase relative to the mean-subtracted HTI YME when the material symmetry is as defined by a set of vertical fractures because torque-on-bit will be loading perpendicular to the axis of material symmetry. In this example the cross-overs are not always synonymous as can be expected for real rock where a continuum of mechanical rock properties will occur based on the natural heterogeneity of a complex, natural system where it is understood that the detection of fractures from among various relationships between these cross-overs is just one implementation among others.

Among the advantages of the method is to use the values of the differences in the elastic coefficients between the various types of cross-over that can be expected from the method where in simple cases fracture flags identified by the simultaneous cross of both YME and PR curves as is afforded by the method, to other cases of mechanical rock heterogeneity that is caused by the occurrence of only one curve crossing over the other or vice versa. Here this may lead to additional fracture classification schemes that would involve the cross-over of one set of coefficients relative to the other.

Elastic coefficients and the variation in the elastic coefficients where the variation of the elastic coefficients is determined in one particular application through the differences as obtained by a subtraction and of the elastic coefficients for the HTI and VTI instances as is specified by the orientation of the well in relation to the axis of material symmetry. When the variations determined from the stress-strain relationships based on the geographical orientation of the well relative to the bedding planes provide an indication of fractures in accordance with the manner described here would provide useful information for the design the emplacement of hydraulic fractures treatments and the selection of hydraulic fracture initiation points.

Other Improvements as are Envisioned by the Method

In one embodiment, the wellbore is drilled laterally through an unconventional shale reservoir. The natural variations in the strength of an unconventional shale reservoir could be viewed by plotting a plethora of stress-strain relationships that could be derived in association with every single turn of the bit (FIG. 8B). This would create a scatter plot that can be analyzed by using statistical methods to find significant relationships in the data that are used to classify the nature of the rock deformation based on the groupings of the mechanical rock property data on the stress-strain diagram.

These classified rock strength measurements may be indexed according to their spatial location along the well bore trajectory. Classification of the type of rock deformation such as strong or weak, brittle or ductile based on the mechanical rock properties and in particular the elastic coefficients in the stress-strain diagrams when logged using the MWD system could be used to identify and select zones along the wellbore selection of hydraulic fracture initiation points for the emplacement of hydraulic fractures.

Other stress-strain relationships can be developed in the manner described here where the forces acting in the formation in connection with the bit are the effective stress obtained by differencing the force and the fluid pressures in the drilling system or the Torque acting on the bit. In another embodiment the orientation and geometry of the cutters in relation to the WOB and Torque can be used to describe tractions that are normal and tangential to the cutting face and can be used to specify additional coefficients of elasticity as provided by the method. The normal traction can be modified by the drilling fluid pressures to calculate effective normal stress acting on the cutting face of the rock formation. When the effective normal stress and shear stress can be projected onto a fault place where the fault plane undergoes reactivation as evidenced by the methods here, when used in conjunction with a failure criterion, can provide critical information on the state-of-stress in the reservoir.

Inclination Measurements as Provided by the ZFL Measurements

Referring again to FIG. 6, it is understood that changes in the inclination of the bit, detectible with MWD data in combinations of the various forms discussed herein, may also identify discontinuities along a bore hole. More specifically, as illustrated, if the drill bit encounters a mechanical discontinuity or geological boundary, the cutting face of the drill bit may change its orientation is response to the orientation and stresses acting on the heterogeneity. Instantaneous changes in inclination of the drill bit relative to a long-term average inclination can be obtained using a sensor or an array of sensors (e.g. axial and lateral or rotary accelerometers) configured to record and extract signals in relation to the three, independent spatial axes of the drilling vibrations. The magnitude of the deflection of the inclination of the drill bit relative to a long-term trend in the direction of the drill bit depends on the orientation of the mechanical discontinuity or geological boundary with respect to the cutting face of the drill bit and therefore the deflection is understood to indicate a change in the mechanical rock properties. If there is an indication of a fracture based on the RMS levels of the other measurements, but no corresponding deflection in the bit, then it is understood that the orientation of the mechanical discontinuity is perpendicular to the orientation of the well bore trajectory. The magnitude of the deflection of the drill bit can be determined through a principle component analysis of the signals extracted from the drilling vibrations where the time window used to obtain principle components of drilling motion of the signals may be normalized by the bit speed and where the principle components can be expressed as changes in the rate of penetration (ROP) or the acceleration of the drill bit.

Discrete Microearthquake Detection Method to Identify a Fault

Referring now to FIG. 15 if the forces acting on the formation in connection with the drill bit and drilling fluid system when conducting drilling operations are sufficient to overcome the failure criteria of a pre-existing fault, then the fault will slip or fail. Reactivation of a fault or pre-existing fracture can be evidenced by extracting a signal from the drilling vibrations that is related to a microseismic event with attendant primary, compressional (P) and secondary, or shear (S) arrivals. In the case where the fault is perpendicular to the trajectory of the wellbore the P-wave arrival is related the particle motion parallel to the axis of the drill string and the S- or transverse wave is the particle motion parallel to the lateral and torsion motion of the drill string. Deviations of the particle motion of the P- and S-waves can be used to determine the orientation of the fault relative to the trajectory of the wellbore.

Reactivation of a fault is expected to create much larger signals are typically expected from the acoustic emissions generated by the fracturing of a rock formation in relation to the cutting action of the bit. By looking at the instantaneous amplitude levels relative to a long-term trend, where the temporal windows used to select the instantaneous amplitudes are related to the bit speed and the long-term window is related to the spatial distribution of the faults in the rock formation, it is possible to identify the location where the bit encountered and crossed a fault.

In the special case of a fault reactivation while conducting drilling operations, the sensors deployed on the bottomhole assembly act like an earthquake seismometer where the geophysical signal processing techniques identify the discrete arrivals of P-waves and S-waves with orthogonal particle motions to detect the presence and reactivation of a pre-existing fault. If the orientation of the fault with respect to the orientation and magnitude of the forces acting on the formation in connection with the bit and drilling fluid system can be determined, then this would enable a method to specify in 3-dimensions a failure criterion of the fault.

Geosteering or Real-Time Applications of the Method

Stick-slip drilling behavior causes the bit speed, typically expressed in revolutions, per minute to increase or decrease according to the distribution of forces acting on the formation in connection with the bit. The variations of the bit speed with respect to the forces used to fracture a rock formation to achieve optimum rates of penetration are typically used to describe the efficiency of the drilling operation, where there is an optimum efficiency that maximizes the rate of penetration with respect to the forces acting on the bit

Techniques that can determine the deformation of a rock formation by describing formation fracturing while conducting drilling provide that when these techniques are enabled in real-time implementations using MWD systems and apparatus, they also can be used to “geosteer” the horizontal well at the bit by maintaining the trajectory of the wellbore in the hydrocarbon bearing zone that experiences fracturing generated by the cutting action of the bit in relation to the forces acting on the formation in connection with the bit and drilling system that has bearing on the mechanical rock properties that will enhance the effectiveness of hydraulic fracture emplacements.

Example embodiments described herein regarding the various control methods may be implemented at least in part in electronic circuitry; in computer hardware executing firmware and/or software instructions; and/or in combinations thereof. Example embodiments also may be implemented using a computer program product (e.g., a computer program tangibly or non-transitorily embodied in a machine-readable medium and including instructions for execution by, or to control the operation of, a data processing apparatus, such as, for example, one or more programmable processors or computers). A computer program may be written in any form of programming language, including compiled or interpreted languages, and may be deployed in any form, including as a stand-alone program or as a subroutine or other unit suitable for use in a computing environment.

FIG. 16 is a flowchart illustrating one method conforming with aspects of the present disclosure. It should be recognized that the elaborate detail and various alternatives and embodiments set out herein may constitute other methods, alone or in combination with that set forth in FIG. 16. Moreover, various steps of the method of FIG. 16, as well as other methods, may be performed within a computer system such as set out in FIG. 17 or may be performed, in whole or in part, in a bottom hole assembly associated with or proximate a bit, such as shown in FIG. 2, may form or be used in steering and may therefore be deployed in the geosteering system illustrated in FIG. 2 or may be deployed in various mechanisms associated with completions. Referring to FIG. 16, the method involves receiving acoustical signals obtained from one or more sensors positioned on a component of a bottom hole assembly (operation 1610). The sensors (e.g., accelerometers) are in operable communication with at least one data memory to store the acoustical signals (e.g., vibration data) where the acoustical signals are generated from a drill bit interacting with a rock formation while drilling a wellbore. The method further involves processing the acoustical signals to obtain at least one set of data values representative of a mechanical rock property of the rock formation along the wellbore created by the drill bit interacting with the rock formation for a period of time (operation 1620). After which, the method may involve identifying a change in the at least one set of data values, where the change is representative of the drill bit crossing a mechanical rock property discontinuity while drilling the wellbore (operation 1630). In some instances, the method may further involve using the data to complete the well (operation 1640).

MRPA and Material Anisotropy or Anisotropy of the Rock Formation

The present disclosure further involves processing drilling dynamics data generated by the deformation of a rock formation in response to forces acting on the rock formation in relation to the alignment of the drill bit. During the drilling process the orientation of the MWD measurements may change with respect to the bedding or material symmetry direction of the rock formation. The constitutive elastic equations or matrix representations of other mechanical properties can be rotated in relation to the orientation of the bit to process the data in a consistent manner with respect to the material symmetry. The processed data may then be used to generate various possible mechanical rock properties.

Referring to FIG. 17, and as discussed extensively above, a method of characterizing rock properties of a formation, based on data gathered while drilling a well bore in the formation and from sensors proximate a drill, first involves receiving and storing signals representative of motions and/or forces on the bit (operation 1710). More particulary, the method may involve receiving acoustical and/or force signals obtained from at one or more sensors positioned on a component of the BHA. The sensors, discussed herein, are in operable communication with at least one data memory to store the signals. The data is initially in the time domain. At the tool or afterward, a processor may transform the acoustical and force signals to obtain spectral estimates of the bit displacements and the forces on the bit (operation 1720). The transformation may involve obtaining an RMS of an amplitude of displacement and amplitude of force.

Many rock formations, such as shale, exhibit anisotropic behavior with respect to their mechanical properties. The rotation of the drill bit and trajectory of the wellbore can change the direction in which the force from the bit is applied to the rock. To align the material symmetry of the rock with the applied force, matrix transformations may be used to rotate the stresses and strains or the material constants, to a common reference frame (operation 1740). Prior to alignment, the spectral estimates may be scaled to obtain the stresses and strains (operation 1730). We describe use cases for rotation of the elastic tensor, but similar principles can be used to rotate other rheological models such as viscoelastic, elastoplastic, visco-elastoplastic, etc.

By processing the stress-strain relationships with respect to the applied force direction, the methods and systems discussed herein can obtain mechancial rock properties of the formation including elastic properties additional to the elastic coefficients discussed above with reference to FIGS. 1-16. For example, in the case of transversely isotropic media, as discussed above, it is possible to obtain the elastic coefficients of C₁₁, C₁₂, C₁₃ and C₃₃; however, the stiffness coefficient C₅₅ may not be obtainable if it is assumed the forces are applied along principle axes and the material symmetry aligns with the principle axes. With rotation, C₅₅ is obtainable.

The five elastic coefficient solution further allows for the computation of the relevant Thomsen parameters used in seismic processing algorithms and provides a representation of the vertical Vp/Vs ratio that has historically been used in vertical wells for reservoir characterization and lithology indication.

General Elastic Tensor Rotation

Rotation of the elastic tensor can be accomplished with Bond transformations. We discuss the rotations in terms of the compliance tensor, but there are many manifestations of Hooke's Law with different elastic coefficients, any of which may be related back to the compliance tensor. The stiffness tensor (C_(ij)'s) is the most common manifestation, and relates to the compliances (S_(ij)'s) as:

S_(ij)=[C_(ij)]⁻¹

The Bond transformation is given by:

S′_(ij)=NS_(ij)N^(T)

where S′_(ij) is the rotated compliance tensor and N is the Bond rotation matrix. In a general Cartesian coordinate system there are three potential Bond transformations (FIG. 18) for rotation about each axis. For rotation about three arbitrary axes 1, 2, and 3 by angles θ, φ, and φ, respectively the Bond transformations for the 6×6 compliance tensor are:

$\quad{\begin{bmatrix} 1 & 0 & 0 & 0 & 0 & 0 \\ 0 & {\cos^{2}\theta} & {\sin^{2}\theta} & {2\cos \; \theta \; \sin \; \theta} & 0 & 0 \\ 0 & {\sin^{2}\theta} & {\cos^{2}\theta} & {{- 2}\cos \; {\theta sin}\; \theta} & 0 & 0 \\ 0 & {{- \cos}\; {\theta sin\theta}} & {\cos \; {\theta sin}\; \theta} & {{\cos^{2}\theta} - {\sin^{2}\theta}} & 0 & 0 \\ 0 & 0 & 0 & 0 & {\cos \; \theta} & {{- \sin}\; \theta} \\ 0 & 0 & 0 & 0 & {\sin \; \theta} & {\cos \; \theta} \end{bmatrix},\begin{bmatrix} {\cos^{2}\varphi} & 0 & {\sin^{2}\varphi} & 0 & {2\cos \; \varphi \; \sin \; \varphi} & 0 \\ 0 & 1 & 0 & 0 & 0 & 0 \\ {\sin^{2}\varphi} & 0 & {\cos^{2}\varphi} & 0 & {{- 2}\cos \; \varphi \; \sin \; \varphi} & 0 \\ 0 & 0 & 0 & {\cos \; \varphi} & 0 & {{- \sin}\; \varphi} \\ {{- \cos}\; \varphi \; \sin \; \varphi} & 0 & {\cos \; {\varphi sin}\; \varphi} & 0 & {{\cos^{2}\varphi} - {\sin^{2}\varphi}} & 0 \\ 0 & 0 & 0 & {\sin \; \varphi} & 0 & {\cos \; \varphi} \end{bmatrix},{{and}\begin{bmatrix} {\cos^{2}\phi} & {\sin^{2}\phi} & 0 & 0 & 0 & {2\cos \; \phi \; \sin \; \phi} \\ {\sin^{2}\phi} & {\cos^{2}\phi} & 0 & 0 & 0 & {{- 2}\cos \; \phi \; \sin \; \phi} \\ 0 & 0 & 1 & 0 & 0 & 0 \\ 0 & 0 & 0 & {\cos \; \phi} & {\sin \; \phi} & 0 \\ 0 & 0 & 0 & {{- \sin}\; \phi} & {\cos \; \phi} & 0 \\ {{- \cos}\; \phi \; \sin \; \phi} & {\cos \; \phi \; \sin \; \phi} & 0 & 0 & 0 & {{\cos^{2}\phi} - {\sin^{2}\phi}} \end{bmatrix}}}$

Any arbitrary rotation can be achieved using a combination of the three Bond transformations.

One use for Bond transformations is to rotate the compliance tensor into a reference frame such that the principle stresses and strains align with the reference frame. Alternatively, the compliance tensor can be left in the reference frame that aligns the material symmetry with the principle axes, and the principle stresses and strains can be rotated to the material symmetry reference frame through other matrix transformations. We only demonstrate rotation of the compliance tensor with Bond transformations in our examples with the understanding that both solutions or any mixture of the two use matrix transformations and produce equivalent results.

Rotation to Account for Wellbore Inclination

In one specific implementation, the compliance tensor is rotated based on the inclination of the wellbore (FIG. 19). The material symmetry is assumed to align with the formation layering and the angle of rotation may be determined by a directional survey combined with the formation dip. The formation dip is determined by a geosteering interpretation in which the gamma ray and possibly other observations while drilling such as mud logs, rate of penetration, etc. are used to infer the position within the reservoir. By mapping out the position within the reservoir along the well, the angle at which the bit is contacting the formation can be determined.

FIG. 19 illustrates a drill bit striking a formation at an oblique angle. The rock material symmetry is indicated by the straight lines for a transversely isotropic model. The vertical direction is aligned with the axis of rotational symmetry, but in general the angle between the axis of rotational symmetry and the drill bit axis represents the angle of rotation. If the well trajectory is set to the 1-direction in FIG. 18, the angle of rotation to account for changes in the inclination of the drill bit would be in the 2-direction, and controlled by the Bond transformation containing rotation angle φ. Let it be understood that any alternative reference frames one may choose to consider, such as setting the well trajectory to the 2-direction, are equivalent forms of the solution we have provided. One may also set the well trajectory to the 1-direction and rotate about the 3-direction. The obtained compliances may be different depending on the choice of reference frame but the predicted strain of the formation to an applied stress would be equivalent.

Rotation to Account for Drill Bit Rotation

In another implementation, the compliance matrix is rotated in accordance with the rotation of the drill bit. Revolutions of the bit can be measured during drilling, for example, with gyros, magnetometers, accelerometers, etc. In the special case of a horizontal wellbore with the wellbore axis aligned in the 1-direction, the bit will apply stress to the formation at every angle in the 2-3 plane during rotation. Assuming the formation follows a transversely isotropic elastic model with the plane of isotropy aligned with the 1-2 plane, rotation of the bit will apply stress at every elastically unique angle and allow for determination of the entire compliance tensor.

One embodiment of the method includes integrating over partial turns of the bit for input stresses and strains that represent finite revolutions of the bit (FIG. 20). For example, when the wellbore is aligned with the 1-axis, partial rotation about the 1-axis can be described in integral form as:

${\int\limits_{\theta_{1}}^{\theta_{2}}{\left( {{{NS}_{ij}N^{T}\sigma_{j}} - e_{i}} \right){\theta}}} = 0$

An analytical expression can be derived for partial rotations by directly substituting the Bond transformation into the integral and evaluating the integral with respect to a starting revolution θ₁ and an ending revolution θ₂. The expression can then be populated with stresses and strains determined from drilling dynamics data and changes in bit revolution to solve for compliances of the rock formation. Here we assume the absolute position of the cutter arms in relation to the material symmetry is known. In a later section we will describe how the absolute position may be determined without independent knowledge of the cutter position.

FIG. 20 illustrates partial rotations of a drill bit, the partial rotations illustrated by θ₁, θ₂, θ₃, and θ₄. Each arrow indicates limits of integration that would be applied. The end of each partial rotation becomes the beginning of the next partial rotation. Integration of the compliances inherently assumes arithmetic averaging of the compliances. Arithmetic averaging of compliances represents the minimum bound on the estimated elastic properties. The maximum bound can be obtained by performing the Bond transformations and integration on stiffnesses or by inverting the compliances and Bond transformations prior to integration. With the two bounding conditions, averaging schemes may be applied to represent the transfer of stress through the rock fabric.

A degenerative case of the partial rotation processing is to integrate over one complete rotation. This simplified case may be useful in areas that have high rotational speed such that partial rotations average over the rock properties are too much to isolate the angular contributions. By averaging over one complete turn of the bit, sensitivity to one or more of the elastic coefficients may be lost.

Another embodiment of the method includes populating the rotated compliance matrix directly with stresses and strains assuming infinitesimal rotations of the bit. As an example, we may separate the drilling dynamics data into incremental rotations of the bit in time, then apply a Fourier transform to each incremental rotation to convert the data into the frequency domain and process the stress-strain relations in a given frequency band. We may want to process the data with this technique if the bit rotation is fast in relation to the Fourier transform window length.

The wellbore inclination and bit rotation can be simultaneously taken into account by first rotating the compliance matrix for the inclination, then applying any embodiment described herein to account for bit rotation.

Arbitrary Rotation to Solve for Material Symmetry and Reference Frame

The previous implementations have assumed a material symmetry and reference frame prior to processing the stress-strain relations. However, in some applications we may not know the material symmetry and reference frame, and we may want to search for these properties during processing. One such instance would be in a fractured reservoir. If the rock fluctuates from an isotropic elastic model to a transversely isotropic model or a transversely isotropic model to orthorhombic or triclinic model, we could search for a representative elastic model and reference frame that best described the stress-strain relations. We always expect a better prediction of the measured data with lower symmetry models due to the increase in parameters. When searching for an elastic model we evaluate the model based on the improvement in our prediction of the measured data. The reference frame has the same number of parameters regardless of the elastic model. We can evaluate the reference frame based on how well that reference frame predicts the measured data.

As an example, if we are implementing the integral solution for partial turns of the bit we may want to dynamically solve for the bit position with respect to the material symmetry. Assuming we know the number of bit revolutions from an auxiliary data set, we can solve for the absolute bit position by:

-   -   a. incrementally shifting the start position (from 0 to 2π)     -   b. populating the stress-strain relations and solving for the         compliance tensor     -   c. using the compliance tensor and stress measurements to         predict the strain measurements (or using the strain         measurements to predict stresses)     -   d. choosing the shift that best predicted the strain         measurements

This approach has been applied to an unconventional shale formation. The well was assumed to be contacting the formation at 90°, or parallel to bedding such that no inclination rotation was required. The anisotropic elastic variations are commonly described in terms of Thomsen coefficients, which are used by seismic processors to represent angular variations in velocity. Here Thomsen coefficients are useful because they have unitless terms that represent the degree of anisotropy. The three unitless Thomsen parameters in terms of stiffness coefficients are:

$ɛ = \frac{C_{11} - C_{33}}{2C_{33}}$ $\gamma = \frac{C_{66} - C_{55}}{2C_{55}}$ $\delta = \frac{\left( {C_{13} + C_{55}} \right)^{2} - \left( {C_{33} - C_{55}} \right)^{2}}{2{C_{33}\left( {C_{33} - C_{55}} \right)}}$

As mentioned previously, if we assume the stresses are applied along the principle axes and the principle axes are aligned with the material symmetry we have no sensitivity to C55. Based on these equations we would only recover one Thomsen parameter. The relative relationships follow the available data from the literature and are within the range of core measured anisotropy parameters.

Stratigraphic Grouping of Similar Rocks to Solve for Elastic Properties, Material Symmetry, or Reference Frame

Another implementation of the method uses nearby wells or the same well passing through a formation multiple times to process stratigraphically similar units of rock. The drilling dynamics data is normally grouped based on the measured depth (MD) or total length encountered while drilling a single well. However, in structurally complex environments adjacent rock may contain highly variable properties due to faulting, or the well may pass through many different formations or faces of rock. The processing may be more robust by organizing the data by stratigraphic section rather than MD. For the special case of wellbore inclination, we may improve the anisotropic property determination by including multiple inclinations for a single rock unit.

If the formation has zero dip and no faults, the well can be sorted according to the true vertical depth (TVD). A dip correction can be applied to the wellbore inclination to convert TVD to true stratigraphic depth (TSD). In the presence of faults, the throw of the fault must be determined to shift the stratigraphy in addition to the dip correction, then the data can be sorted by TSD. To process partial rotations in TSD the data is sorted by TSD, data for processing is selected, then the data is resorted back to MD for sequential processing of each section. For partial rotations the input data must be sequential, but for inclination processing only, the sequence can be mixed.

FIG. 21 is a block diagram of a machine in the example form of a computer system 2100 within which instructions 2106 for causing the machine to perform any one or more of the methodologies discussed herein may be executed by one or more hardware processors 2102. In various embodiments, the machine operates as a standalone device or may be connected (e.g., networked) to other machines. In a networked deployment, the machine may operate in the capacity of a server or a client machine in server-client network environment, or as a peer machine in a peer-to-peer (or distributed) network environment. Further, while only a single machine is illustrated, the term “machine” shall also be taken to include any collection of machines or controllers that individually or jointly execute a set (or multiple sets) of instructions 2106 to perform any one or more of the methodologies discussed herein, including that set out in FIG. 16 as well as the various methodologies discussed herein to obtain and/or compute ZFL, axial displacement, rotary displacement, axial and rotary accelerations, bit displacement, RMS values of various measurements, Poisson's ratio, Young's modulus of elasticity, and others.

As depicted in FIG. 21, the example computing system 2100 may include one or more hardware processors 2102, one or more data storage devices 2104, one or more memory devices 2108, and/or one or more input/output devices 2110. Each of these components may include one or more integrated circuits (ICs) (including, but not limited to, field-programmable gate arrays (FPGAs), application-specific ICs (ASICs), and so on), as well as more discrete components, such as transistors, resistors, capacitors, inductors, transformers, and the like. Various ones of these components may communicate with one another by way of one or more communication buses, point-to-point communication paths, or other communication means not explicitly depicted in FIG. 21. Additionally, other devices or components, such as, for example, various peripheral controllers (e.g., an input/output controller, a memory controller, a data storage device controller, a graphics processing unit (GPU), and so on), a power supply, one or more ventilation fans, and an enclosure for encompassing the various components, may be included in the example computing system 2100, but are not explicitly depicted in FIG. 21 or discussed further herein.

The at least one hardware processor 2102 may include, for example, a central processing unit (CPU), a microprocessor, a microcontroller, and/or a digital signal processor (DSP). Further, one or more hardware processors 2102 may include one or more execution cores capable of executing instructions and performing operations in parallel with each other. In some instances, the hardware processor is within the bit sub, and others it is part of another separate processing system.

The one or more data storage devices 2104 may include any non-volatile data storage device capable of storing the executable instructions 2106 and/or other data generated or employed within the example computing system 2100. In some examples, the one or more data storage devices 2104 may also include an operating system (OS) that manages the various components of the example computing system 2100 and through which application programs or other software may be executed. Thus, in some embodiments, the executable instructions 2106 may include instructions of both application programs and the operating system. Examples of the data storage devices 2104 may include, but are not limited to, magnetic disk drives, optical disk drives, solid state drives (SSDs), flash drives, and so on, and may include either or both removable data storage media (e.g., Compact Disc Read-Only Memory (CD-ROM), Digital Versatile Disc Read-Only Memory (DVD-ROM), magneto-optical disks, flash drives, and so on) and non-removable data storage media (e.g., internal magnetic hard disks, SSDs, and so on).

The one or more memory devices 2108 may include, in some examples, both volatile memory (such as, for example, dynamic random access memory (DRAM), static random access memory (SRAM), and so on), and non-volatile memory (e.g., read-only memory (ROM), flash memory, and the like). In one embodiment, a ROM may be utilized to store a basic input/output system (BIOS) to facilitate communication between an operating system and the various components of the example computing system 2100. In some examples, DRAM and/or other rewritable memory devices may be employed to store portions of the executable instructions 2106, as well as data accessed via the executable instructions 2106, at least on a temporary basis. In some examples, one or more of the memory devices 2108 may be located within the same integrated circuits as the one or more hardware processors 2102 to facilitate more rapid access to the executable instructions 2106 and/or data stored therein.

The one or more data storage devices 2104 and/or the one or more memory devices 2108 may be referred to as one or more machine-readable media, which may include a single medium or multiple media that store the one or more executable instructions 2106 or data structures. The term “machine-readable medium” shall also be taken to include any tangible medium that is capable of storing, encoding, or carrying instructions 2106 for execution by the machine and that cause the machine to perform any one or more of the methodologies of the present invention, or that is capable of storing, encoding, or carrying data structures utilized by or associated with such instructions 2106.

The input/output devices 2110 may include one or more communication interface devices 2112, human input devices 2114, human output devices 2116, and environment transducer devices 2118. The one or more communication interface devices 2112 may be configured to transmit and/or receive information between the example computing system 2100 and other machines or devices by way of one or more wired or wireless communication networks or connections. The information may include data that is provided as input to, or generated as output from, the example computing device 2100, and/or may include at least a portion of the executable instructions 2106. Examples of such networks or connections may include, but are not limited to, Universal Serial Bus (USB), Ethernet, Wi-Fi®, Bluetooth®, Near Field Communication (NFC), and so on. One or more such communication interface devices 2112 may be utilized to communicate one or more other machines, either directly over a point-to-point communication path or over another communication means. Further, one or more wireless communication interface devices 2112, as well as one or more environment transducer devices 2118 described below, may employ an antenna for electromagnetic signal transmission and/or reception. In some examples, an antenna may be employed to receive Global Positioning System (GPS) data to facilitate determination of a location of the machine or another device.

In some embodiments, the one or more human input devices 2114 may convert a human-generated signal, such as, for example, human voice, physical movement, physical touch or pressure, and the like, into electrical signals as input data for the example computing system 2100. The human input devices 2114 may include, for example, a keyboard, a mouse, a joystick, a camera, a microphone, a touch-sensitive display screen (“touchscreen”), a positional sensor, an orientation sensor, a gravitational sensor, an inertial sensor, an accelerometer, and/or the like.

The human output devices 2116 may convert electrical signals into signals that may be sensed as output by a human, such as sound, light, and/or touch. The human output devices 2116 may include, for example, a display monitor or touchscreen, a speaker, a tactile and/or haptic output device, and/or so on.

The one or more environment transducer devices 2118 may include a device that converts one form of energy or signal into another, such as from an electrical signal generated within the example computing system 2100 to another type of signal, and/or vice-versa. Further, the transducers 2118 may be incorporated within the computing system 2100, as illustrated in FIG. 21, or may be coupled thereto in a wired or wireless manner. In some embodiments, one or more environment transducer devices 2118 may sense characteristics or aspects of an environment local to or remote from the example computing device 2100, such as, for example, light, sound, temperature, pressure, magnetic field, electric field, chemical properties, physical movement, orientation, acceleration, gravity, and so on. Further, in some embodiments, one or more environment transducer devices 2118 may generate signals to impose some effect on the environment either local to or remote from the example computing device 2100, such as, for example, physical movement of some object (e.g., a mechanical actuator), receiving or processing accelerometer data, strain gauge data, and the like.

SUMMARY

As such, under the aforementioned considerations the spatial variations in the MWD data in in particular the MWD vibrations as obtained through a combination of one or more of the measurements are used to describe variations in mechanical rock properties where the variations and occurrences of the can be described as fractures based on the methods used here. Here calculations made using signal processing techniques are used to populate innovative new stress-strain relationships.

The method discloses how the spatial variations in one or more or a combination of the mechanical rock properties as can be obtained from the stress-strain relationships used to identify the nature and occurrence of fractures, fracture swarms and other mechanical discontinuities (boundaries) such as bedding planes and/or faults that offset or otherwise separate rock formations with different mechanical rock properties.

The present disclosure uses an innovative, new methodologies to determine the deformation of a rock formation by systematically relating forces acting on a rock formation in connection with the drill bit and drilling fluid system to the geophysical signal processing of drilling vibrations generated by the fracturing of the rock in response to the cutting action of the bit to obtain a strain measurement. This approach allows elastic coefficients to be derived that can be used, in general, to describe the general deformation of a rock formation in response to the forces acting on a rock formation.

The elastic coefficients can be used to determine the deformation of a hydrocarbon bearing formation in response to the forces acting on the formation during the emplacement of hydraulic fractures in connection with hydraulic fracture stimulation treatments and in particular based on the spatial variations of the elastic coefficients select hydraulic fracture initiation points for the purpose of enhancing the effectiveness of the hydraulic fracture emplacements in low permeability hydrocarbon bearing rock formations.

Methods that can in general understand and predict the nature and occurrence of fracturing along a wellbore are important considerations for the specification and design of hydraulic fracture stimulation techniques because decisions to either eliminate zones that do not experience brittle deformation or select zones that exhibit more brittle deformation relative to other zones can be used to improve the economics of the production. Further where the rock deformation is accommodated by failure of pre-existing fractures, these pre-existing fractures should be avoided when setting packers during the pumping operations or they can be targeted for hydraulic fracture stimulation depending on the desired approach to be used when completing the well.

The following is a list, albeit not exhaustive or limited to only the list, of innovations.

A method of processing drilling vibration information to obtain one or more measurements such as the RMS acceleration or the sizes and displacements of the fractures where the measurements are understood to correspond to the motion of the bit and in particular where the motion of the bit is expressed in the ZFL of the displacement spectra to inform the depth of cut as a displacement per revolution or displacement per time to make a revolution. Therefore, spatial variation in the of the motion of the bit in relation to the forces acting on the bit are systematically investigated using stress-strain constitutive equations where the variables are populated using geophysical signal processing of the MWD data to obtain the forces acting on the bit and the motion of the bit in relation to the forces in order to solve for the elastic coefficients where the elastic coefficients are the YME and PR and where the elastic coefficients are determined using stress-strain constitutive relationships that represent the case of where the drilling is perpendicular to an axis of material symmetry and parallel to an axis of material symmetry when the elastic coefficients are taken to represent spatial variations in mechanical rock properties where the spatial variations are used in predictable ways such are used to identify fractures, fracture swarms and other mechanical discontinuities (boundaries) such as faults and bedding planes that offset or otherwise act to separate rock formations with different rock properties.

A further method of processing drilling vibrations in relation to the forces acting on a rock formation in connection with a drill bit and drilling fluid system to specify stress-strain relationships to derive elastic coefficients that in general predict the deformation of a rock formation and in particular, in one form are used to predict the deformation of a hydrocarbon bearing formation in response to the forces acting on the hydrocarbon bearing formation where the forces are generated by the emplacement of hydraulic fractures along a horizontal well to enhance the effectiveness of hydraulic stimulation treatments in a hydrocarbon bearing rock formations generally involves the following:

-   -   (1) Deploying sensors or an array of sensors that are sensitive         to the amplitudes and frequencies of the drilling vibrations         generated in response to the fracturing in relation to the         cutting action of the bit using an instrumented sonde attached         to a bottomhole assembly to record one or more axis of the         drilling vibrations and where the sonde may employ sensors to         record the drilling efficiencies such as the bit speed and the         forces acting on the formation in connection with the drill bit         and fluid pressure changes created by the circulation of the         drilling fluids while conducting drilling operations.     -   (2) Extracting signals from the drilling vibrations that are         related to the elastic energy released by deformation and         failure of the rock formation in response to the cutting action         of the bit that are transmitted from the rock formation in         connection with the bit to the receivers through the drilling         collar to the instrumented sonde in the bottomhole assembly by         filtering the unwanted noise and other interfering modes of wave         propagation not related to the tool arrival, generally involves         using at least one or a combination of filters through the         application and use of (i) mechanically isolated sensors from         energy transmitted through the rock formation and drilling fluid         system (ii) bandpass or FK spatial filters, (iii) multiple         receivers to stack the signals and attenuate the random noise         or (iv) placing the sonde in a position in the drill string that         is sufficiently far enough away from the bit to attenuate the         higher frequencies of the unwanted wave arrivals.     -   (3) Using geophysical signal processing techniques to measure         the temporal and spatial variations in the amplitude and         frequency of the signals as can be extracted or otherwise         derived from the drilling vibrations to determine the nature and         occurrence of fracturing, fracture swarms and other mechanical         discontinuities (boundaries) such as faults and bedding planes         that offset or otherwise act to separate rock formations with         different rock properties using one or a combination of:         -   a. RMS measurements of the signals extracted from the             drilling vibrations where the signals may be expressed in             terms of displacement, vibration or acceleration, and where             the time windows used for the RMS calculation may be             normalized by the bit speed in order to use the relative             changes in the RMS level to identify locations where the             changes in the instantaneous RMS levels or the RMS level of             the bit displacement, velocity or acceleration in connection             with at least one revolution of the bit relative to the RMS             level averaged over longer spatial and temporal interval are             described to identify in the mechanical rock properties as             are expected in relation to the bit encountering or crossing             a mechanical boundary or other geological discontinuity.         -   b. Using geophysical signal processing techniques to measure             the temporal and spatial variations in the amplitude and             frequency on at least one component of the signals extracted             from the drilling vibrations using methods that are             generally recognized as appropriate for the analysis of             microearthquake source mechanisms in order to measure the             size and displacements of the fracturing in relation to the             cutting action of the bit, where the variations in the size             and displacements of the fracturing are used to indicate the             presence of pre-existing faults or fractures or the presence             of other mechanical or geological boundaries as encountered             by the drill bit.         -   c. Instantaneous change in inclination when using a sensor             to record and extract signals related to the three,             independent spatial axes of the drilling vibrations to             identify the fracture or bedding plane orientation relative             to the orientation of the drilling face, where the drilling             face is orthogonal to the inclination of the principle axial             vibration as can be viewed through a principle component             analysis of signals extracted from the drilling vibrations             where the time window used to obtain average of the signals             may be normalized by the bit speed.         -   d. Instantaneous change in forces acting on the bit when             using a sensor to record and extract signals related to the             forces acting on the bit where the time window used to             obtain average of the signals may be normalized by the bit             speed.         -   e. Reactivation of a fault or pre-existing fracture as             evidenced by extracting a signal from one or more of the             principle components of the drilling vibrations that is             related to a microseismic event with attendant primary,             compressional (P) and secondary, or shear (S) arrivals,             where the P-wave arrival is related the particle motion             parallel to the axis of the drill string and the S- or             transverse wave is the particle motion parallel to the             lateral and torsion motion of the drill string.         -   f. Deriving empirical relationships that relate the             variations in the RMS measurements to describe the nature             and occurrence of fractures by establishing multivariate             relationships between the aforementioned measurements with             fracture imaging logs, and mud logs and other petrophysical             logs where the measurements as provided by the logs are             sensitive to the presence of a mechanical boundary or other             geological discontinuity.         -   g. In another embodiment, the signals are used to derive             empirical relationships that relate the variations in the             measurement of one or more of the signal with fracture             imaging logs, and mud logs and other petrophysical logs or             instantaneous drilling dynamics such as ROP where the             measurements as provided by the logs are sensitive to the             presence of a mechanical boundary or other geological             discontinuity. These relationships provide a method to             classify the nature and occurrence of fractures by             establishing multivariate relationships between the             aforementioned measurements and in particular where the             measurements involve the differences between elastic             coefficients as are determined by the stress-strain             relationships in relation to the orientation of drilling             with respect to the axis of material symmetry.     -   (4) A method using geophysical signal processing techniques to         systematically relate measurements of the forces acting on rock         formation in connection with the drill bit and drilling fluid         system (stress) to the variations in the fracturing of a rock         formation in response to the cutting action of the bit (strain)         to obtain innovative, new stress-strain relationships where the         application and use of the stress-strain relationships allow for         the derivation of elastic coefficients where the stress-strain         relationships involve a combination of stress and strain from         one or more of the following:         -   a. The strain is the zero frequency level of bit             displacement, be it axial or and where the axial or lateral             component of motion may be determined from a principle             component analysis of the bit motion and where the time             window used to measure the RMS level is proportional to the             bit speed.         -   b. The strain is the displacement on the fractures generated             in response to the cutting action of the bit provided by             signal processing techniques in relation to the application             and use of models used to describe microearthquake source             parameters.         -   c. The strain is a volumetric strain determined through a             combination of the displacement of the fracturing and the             area of the fracturing in response to the cutting action of             the bit as is provided by signal processing techniques in             relation to the application and use of models used to             describe microearthquake source parameters.         -   d. Where the stress is determined from measurements taken             while drilling in relation to (i) the forces acting on the             formation in connection in connection with the bit and fluid             system, such as the weight on the bit (WOB) or torque-on-bit             or (ii) the rms amplitude of the acceleration of the bit.     -   (5) A method to generally derive elastic coefficients using         stress-strain relationships determined by the cutting action of         the bit in relation to the forces acting on the formation in         connection with the bit and drilling fluid system to:         -   a. predict the deformation of a hydrocarbon bearing             formation in response to the forces acting on the             hydrocarbon bearing formation where the forces are generated             by the emplacement of hydraulic fractures in connection with             a hydraulic fracture stimulation treatment along a             horizontal well.         -   b. assist, in real-time, the steering of a drilling bottom             hole assembly (BHA) in order to maintain the tracking of the             drill bit through geological formations as are targeted             according to the desired mechanical rock properties as             specified through the elastic coefficients, especially where             the mechanical rock properties are relevant to the             production of commercially significant hydrocarbons using             hydraulic fracture stimulation techniques.         -   c. Classify the elastic coefficients on the stress-strain             diagrams to identify zones of brittle rock to target             fracture initiation points         -   d. Geologically correlate the mechanical rock properties of             a formation using the elastic coefficients or a use a             classification of the elastic coefficients such as brittle             and ductile or strong and weak between two to predict the             nature and occurrence of mechanical rock properties in the             locations of undrilled wells         -   e. Using the variation of the elastic coefficients in             predictable ways to identify locations in the mechanical             rock properties where the variations in the mechanical rock             properties are localized and discrete in order to identify             the nature and occurrence of fractures, fracture swarms and             other mechanical discontinuities (boundaries) such as             bedding planes and/or faults that offset or otherwise             separate rock formations with different mechanical rock             properties to:             -   i. Geologically correlate the mechanical discontinuities                 or geological boundaries across two more wells where in                 the case an undrilled well lies within the region that                 is contained by the correlation of two or more wells,                 predict the nature and occurrence of the mechanical                 discontinuities and geological boundaries in an                 undrilled well             -   ii. Locate swell packers away from mechanical                 discontinuities such as fracture swarms or geological                 boundaries to improve pressure isolation in a                 hydrocarbon bearing formation             -   iii. Target natural fractures for initiation of                 hydraulic fractures in connection with the stimulation                 and treatment of a well in a hydrocarbon bearing                 formation             -   iv. Isolate locations in the hydrocarbon bearing                 formation based on the nature and occurrence of                 mechanical discontinuities or geological boundaries to                 induce more fractures in competent reservoir, thereby                 increasing contacted reservoir             -   v. Provide correlation information for seismic                 inversions         -   f. where the wells are drilled in two mutually perpendicular             directions that where the specification of the elastic             coefficients is made using the orientation of the well in             respect to the principle axes of tectonic stress or the             state-of-stress acting on a rock formation         -   g. where the undrilled well is later drilled using             geosteering         -   h. where the techniques used to determine the nature and             occurrence of fracturing, fracture swarms and other             mechanical discontinuities (boundaries) such as faults and             bedding planes that offset or otherwise act to separate rock             formations with different rock properties using one or a             combination of the measurements as provided by the             geophysical signal processing techniques and the elastic             coefficients as provided by claim 4.         -   i. where the drilling efficiencies are varied in a             deliberate and systematic manner to generate an in-situ             stress-strain relationship.

While the present disclosure has been described with reference to various implementations, it will be understood that these implementations are illustrative and that the scope of the disclosure is not limited to them. Many variations, modifications, additions, and improvements are possible. More generally, implementations in accordance with the present disclosure have been described in the context of particular implementations. Functionality may be separated or combined in blocks differently in various embodiments of the disclosure or described with different terminology. These and other variations, modifications, additions, and improvements may fall within the scope of the disclosure as defined in the claims that follow. 

What is claimed is:
 1. A method of characterizing rock properties comprising: obtaining, using a processor, dynamic signals from one or more sensors positioned on a component of a bottom hole assembly, the dynamic signals generated from a drill bit interacting with a rock formation while drilling a wellbore; aligning, using the processor, the dynamic signals in relation to a position of the drill bit with respect to an axis of symmetry of the rock formation; and processing, using the processor, the aligned dynamic signals to obtain at least one set of data values representative of a mechanical rock property of the rock formation along the wellbore created by the drill bit interacting with the rock formation for a period of time.
 2. The method of claim 1 wherein the dynamic signals include acoustic signals generated from vibrations of the drill bit interacting with the rock formation while drilling the wellbore and force signals based on at least one of a torque-on-bit or a weight-on-bit of the drill bit interacting with the rock formation.
 3. The method of claim 1 wherein aligning the dynamic signals in relation to the position of the drill bit with respect to the axis of symmetry of the rock formation includes: transforming the dynamic signals to frequency domain signals; and rotating the frequency domain signals in relation to a position of the drill bit with respect to an axis of symmetry of the rock formation.
 4. The method of claim 3 wherein rotating the dynamic signals in relation to the position of the drill bit with respect to the axis of symmetry of the rock formation includes applying a matrix transformation to rotate a reference frame of the rock formation.
 5. The method of claim 3 wherein rotating the dynamic signals in relation to the position of the drill bit with respect to the axis of symmetry of the rock formation includes applying a matrix transformation to rotate the frequency domain signals.
 6. The method of claim 3 wherein processing the aligned dynamic signals includes processing the aligned frequency domain signals to obtain at least one set of data values representative of a mechanical rock property of the rock formation along the wellbore created by the drill bit interacting with the rock formation for a period of time.
 7. The method of claim 6 wherein the mechanical rock property further includes an elastic coefficient C₅₅.
 8. The method of claim 7 wherein the mechanical rock property includes a plurality of stiffness coefficients of C₁₁, C₁₂, C₁₃ and C₃₃.
 9. The method of claim 1 wherein the mechanical rock property includes a stiffness coefficient C₅₅.
 10. The method of claim 1 further comprising identifying, using the processor, a change in the at least one set of data values, the change representative of the drill bit crossing a mechanical rock property discontinuity while drilling the wellbore.
 11. The method of claim 1 wherein aligning, using the processor, the dynamic signals in relation to the position of the drill bit with respect to the axis of symmetry includes applying a matrix transformation to rotate the dynamic signals, including acoustic or mechanical signals, by applying a Bond transformation.
 12. The method of claim 11 wherein processing the rotated acoustic or mechanical signals comprises obtaining an inclination of the drill bit and rotating a material symmetry axis of the rock formation to align with a trajectory of the drill bit.
 13. The method of claim 11 wherein processing the rotated acoustic or mechanical signals further comprises obtaining a bit rotation speed and rotating the material symmetry axis to account for partial rotations of the bit through integration of the Bond transformation.
 14. The method of claim 11 wherein processing the rotated acoustic or mechanical signals further comprises arbitrarily rotating the material symmetry axis to determine material symmetry orientation.
 15. The method of claim 11 wherein processing the rotated acoustic or mechanical signals further comprises arbitrarily rotating any number of material symmetries to invert for an unknown material symmetry and orientation.
 16. The method of claim 11 wherein processing the rotated acoustic or mechanical signals further comprises generating an estimation of the elastic anisotropy of the material.
 17. The method of claim 11 wherein processing the rotated acoustic or mechanical signals further comprises determining the presence of fractures and/or fracture orientation based on a change in material symmetry.
 18. A method of characterizing rock properties while drilling comprising: receiving accoustical signals obtained from one or more sensors positioned on a component of a bottom hole assembly, the acoustical signals generated from a drill bit interacting with a rock formation while drilling a wellbore; aligning the acoustical signals in relation to a position of the drill bit with respect to an axis of symmetry of the rock formation; processing the aligned acoustical signals to obtain forces acting on the drill bit interacting with the rock formation while drilling the wellbore and to obtain displacements of the drill bit interacting with the rock formation while drilling the wellbore; and processing the forces and the displacements to obtain at least one set of data values representative of a mechanical rock property of the rock formation along the wellbore created by the drill bit interacting with the rock formation for a period of time.
 19. The method of characterizing rock properties of claim 18 further comprises: transforming the motion and force signals to obtain spectral estimates of a displacement of the drill bit and a force on the drill bit; scaling the spectral estimates to obtain stresses and strains; and wherein aligning the motion and force signals in relation to a position of the drill bit with respect to an axis of symmetry of the rock formation comprises rotating a reference frame such the stresses and strains align along a principle axis of the reference frame.
 20. The method of claim 18 further comprising: using a stress strain relationship for processing the forces and displacements to obtain the at least one set of data values representative of a mechanical rock property of the rock formation along the wellbore created by the drill bit interacting with the rock formation for a period of time
 21. The method of claim 18 further comprising: identifying a change in the at least one set of data values, the change representative of the drill bit crossing a mechanical rock property discontinuity while drilling the wellbore.
 22. The method of claim 18 wherein the at least one set of data values includes stiffness coefficients (Cijs).
 23. The method of claim 22 wherein the stiffness coefficients comprise C₁₁, C₁₂, C₁₃, C₃₃, and C₅₅.
 24. The method of claim 18 wherein the one or more sensors are in operable communication with at least one data memory to store the acoustical signals.
 25. The method of claim 18 wherein the acoustical signals are from vibrations generated from the drill bit interacting with the rock formation while drilling the wellbore.
 26. The method of claim 18 wherein the acoustical signals include an axial acceleration of the drill bit and a lateral or rotary acceleration of the drill bit.
 27. The method of claim 18 further comprising: generating a well log spatially identifying mechanical rock property discontinuities along the wellbore using the mechanical rock properties, the mechanical rock property discontinuities indicative of preexisting fractures.
 29. An apparatus comprising: a processing unit in communication with at least one tangible machine readable media including computer executable instructions to perform the operations of: obtaining dynamic signals from one or more sensors positioned on a component of a bottom hole assembly, the dynamic signals generated from a drill bit interacting with a rock formation while drilling a wellbore; aligning the dynamic signals in relation to a position of the drill bit with respect to an axis of symmetry of the rock formation; processing the aligned dynamic signals to obtain at least one set of data values representative of a mechanical rock property of the rock formation along the wellbore created by the drill bit interacting with the rock formation for a period of time.
 30. The apparatus of claim 28, wherein the at least one set of data values includes the shear modulus C₅₅. 